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ANNIVERSARY OF THE SCIENTIST 





Zakovorotny V. L., engineer, scientist, is 80 


12 February, Vilor Lavrentievich Zakovorotny, Honoured Worker of National Sci- 
ence, Dr.Sci. (Engineering), professor, turned 80. 

Vilor Lavrentievich Zakovorotny is a well-known international expert in the field of 
nonlinear dynamics of controlled systems. He has studied the dynamic monitoring and con- 
trol of mechatronic systems and metalworking processes, developed the research framework 
of this area, and created a school of thought. His R&D results were introduced into produc- 
tion at 17 enterprises in Russia and neighboring countries. 

For the past five years, the solutions proposed by V. L. Zakovorotny to the dynamic 
monitoring of machining processes on metal-cutting machines and to the dynamic 
diagnostics of rotor systems have been introduced into the manufacturing practice of 
Rostvertol, Russian Helicopters JSC, Rostselmash Combine Plant LLC, Azov 
Optomechanical Plant JSC, All-Russian Research Institute “Gradient” JSC, and others. 

V. L. Zakovorotny has made a significant contribution to the development of engineering education in Russia, 
improving the educational process, training top-qualification technical and scientific manpower. He has trained more 
than 1,500 engineers, 7 doctors and 42 candidates of science, who work in universities and enterprises of the Rostov 
region. 





V.L. Zakovorotny was one of the originators of training engineering personnel majoring in “Technological- 
Process Automation & Process Control” and “Management and IT in Engineering Systems”. He carries out big job at 
the Academic Methodological Association on Automated Engineering. 

In 1991-2007, as Pro-rector for Research, Prof. Zakovorotny made a great contribution to the reorganization of 
an industry-based university, Rostov Institute of Agricultural Engineering (RIAE) into a comprehensive technical uni- 
versity, DSTU. 

V. L. Zakovorotny is the deputy chairman of the dissertation council D.212.058.02, DSTU, a member of the 
joint dissertation council, SFedU, a member of the editorial boards of three scientific journals included in the list of the 
State Commission for Academic Degrees and Titles of the Russian Federation. Prof. Zakovorotny is the author and co- 
author of 356 publications including seven monographs, five textbooks and teaching aids, 49 patents and copyright cer- 
tificates. 

For his outstanding merits, V. L. Zakovorotny is awarded with the second-class medal of the Order of Merit for 
the Motherland, the titles of Laureate of the State Prize of Ukraine in science and engineering, Honored Scientist of the 
Russian Federation, Honorary Professor and Professor of the year, DSTU. The Russian Air Force Commander awarded 
V. L. Zakovorotny with a Memorial Sign for his contribution to the construction and development of the Russian Air 
Force. 
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Research and modeling of a multilayer composite material using basalt fabric 


A. Karnoub’, D. B. Nezhizhimov’, K. S. Shirinyan* 
"ETH (Zurich, Switzerland) 
*-3Don State Technical University (Rostov-on-Don, Russian Federation) 





Introduction. The range of use of composite materials (CM) is constantly expanding, finding application in many areas 
of mechanical engineering, agricultural technology, aircraft manufacturing, instrumentation, shipbuilding, in the manu- 
facture of high-pressure containers, etc. Quite often, multilayer composites consisting mainly of one type of reinforcing 
material and a binder are used. Of particular interest is the use of various types of reinforcing materials — more durable 
in the places of maximum stress in the cross section — in a single composite. As an example, we can use glass and basalt 
fabrics and fibers using one type of binder. The work objective is to study properties of such a material and to model it 
using the finite element method. 

Materials and Methods. The components used are commercially available. BT-11 basalt fabric, TR-0.5 fiberglass, as 
well as glass mat with a density of 300 g/m’ were used as reinforcing materials. An epoxy resin of the ED-20 grade with 
a PEPA hardener was used as a binder. Two types of material were also manufactured for tensile and bending tests, 
respectively, the differences of which consisted in the number, type and layer sequence. For modeling, CAD COM- 
PASS 3D, APM-FEM module was used. 

Results. Basalt fabric is used in the outer layers of the composite material, fiberglass — in the inner layers. This approach 
provides increasing the tensile strength of the composite during tensile and bending; however, critical failure leads to an 
instant loss of the bearing capacity of the material. The use of glass mat as the core of the composite material showed 
lower allowable stresses, both tensile and bending; but in case of bending, it turned out that when the material was de- 
laminated, the load-bearing capacity of the material was about 10% of the maximum. Modeling of the material is possi- 
ble with some assumptions, in view of the size of the final elements. 

Discussion and Conclusions. The use of basalt fabrics as a reinforcing material provides obtaining products with the 
properties of both glass and carbon plastics. Such a CM will be slightly more expensive than fiberglass and much 
cheaper than carbon fiber. Products made of composite materials (equivalent to isotropic materials) can be modeled in 
computer-aided design systems using the finite element method. It is important to consider the type of loading on the 
product, since CM mainly have anisotropic properties (the load is applied taking into account the direction of fibers). In 
multilayer CM from structural fabrics, it is necessary to direct the loads along the fibers. In addition, it is necessary to 
consider the interlayer shear, different adhesion between the layers, etc. The main assumption of this method is the 


“constancy” of the material thickness, the number of layers and the order of their location. 
Keywords: composite material, basalt fabric, fiberglass, glass mat, finite element method. 
For citation: A. Karnoub, D. B. Nezhizhimov, K. S. Shirinyan. Research and modeling of a multilayer composite material 


using basalt fabric. Vestnik of DSTU, 2020, vol. 20, no. 1, pp. 5-14. https://doi.org/10.23947/1992-5980-2020-20-1-5-14 


Introduction. In the work, polymer composite materials (CM) with a reinforcing composition based on basalt 


and glass fabrics were studied and modeled. Epoxy resin was used as a binder. The 11- and 13-layer composites were 


investigated. 
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In [1], technological and other characteristics of the CM are noted, including those determining the load distri- 
bution in the layers. Here, the principle of “simple” modeling of composite material in various computer-aided design 
(CAD) systems is considered, as well as the behavior of CM depending on the type of loading. 

A method is described for increasing the interlayer shear resistance in polymer CM by adding finely divided 
solid particles (for example, glass) to the binder’. The process of occurrence of interlayer stresses was studied [2]. It is 
possible that the application of this method would avoid stratification and increase the allowable stress of the samples. 

In [3-5], the volume content of the binder, the methods of its application to the reinforcing material, the laying- 
up method, the operating procedure, etc., are determined. 

It is known that the type of deformation and fracture of a composite polymer is determined by its shape, quality 
and binder (including its volume content)’. The corresponding dependence has been derived [6, 7]. 

In [8-10], modeling a part from CM using the finite element method (FEM) in CAD is described. 

The paper presented also uses FEM. With its help: 

- reliability of determining the stresses in the places of destruction of the samples and safety factors is validat- 
ed, 

- specifications of using assumptions that facilitate the calculation of products from polymer CM tested previ- 
ously on tensile machines, are indicated. 

In [11], the influence of the shape and type of reinforcing components on the thermal and mechanical proper- 
ties of a polymer composite material is described. 

In [12, 13], the applicability of composites in engineering, automobile production, construction, etc., were not- 
ed. 

Research Objectives: 

- to justify the use of basaltic fabrics in polymer CM, 

- to determine their strength properties with various reinforcing components, 

- to check the possibilities of modeling CM by the finite element method. 

Materials and Methods 

A composite material based on basalt and glass cloth and glass mat was studied in this work. As a binder, 
epoxy resin ED-20 with hardener polyethylene polyamine (PEPA) was used. The reinforcing agents were fiberglass 
TR-0.5, basalt fabric BT-11 and fiberglass with a density of 300 g/m’. A feature of the experiment was the use of sever- 
al types of reinforcing material in one composite at once. 

Justification for the application of basalt fabric and finite element modeling. At present, CM rarely uses 
basalt fabric, although it possesses characteristics important for high-quality composites: high impact strength, high 
specific strength, heat resistance, environmental cleanliness, high resistance to corrosion and acid, low thermal conduc- 
tivity, affordable price, radio transparency, and good sound absorbing properties. It is also worth noting that basalt is a 
rock, and its reserves are practically unlimited. In the production of basalt fiber, the rock is melted, threads are drawn 
from it, which are used to make fabric, roving, etc.). In terms of strength, basalt fiber is superior to fiberglass and ap- 
proaches carbon fiber. The price of basalt fiber is slightly higher than fiberglass, but significantly lower than carbon 
fiber. This allows fabricating products of higher quality than fiberglass, but at the same time cheaper than carbon fiber. 
Of particular interest is the use of several reinforcing materials in one CM: glass and basalt fabric, as well as glass mat. 
When modeling composites in CAD, it is required to set materials with various characteristics, which creates known 


difficulties. In the framework of this work, the possibility of a standard study of structures by the FEM is determined. 


' Nezhizhimov DB. Sposob uvelicheniya soprotivleniya mezhsloinomu sdvigu v mnogosloinykh kompozitnykh materialakh [Method for increasing 
interlayer shear resistance in multilayer composite materials]. In: Proc.3rd All-Russ. Sci. Conf., Kursk, 2019. P. 84-88. 

* Antibas IR, D'yachenko AG. Issledovanie protsessa razrusheniya sloistogo kompozitnogo materiala [Study on destruction of the layered composite 
material]. In: Proc. 10th Int. Sci.-Pract. Conf. within framework of the 20th Int. Agroindustrial Exhibition “Interagromash-2017”, Rostov-on-Don, 
2017. P. 179-181. 
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For this, anisotropic material is “equated” to isotropic, and the products are designed so that loads can be applied along 
the direction of the fibers. 

The use of structural fabrics makes it possible to model the behavior of CM with some assumption. However, 
it is necessary to consider the type of weaving of fabrics and the angle of rotation between the layers. Thus, it is re- 
quired to obtain the results of an experiment on determining the strength properties, and to find out the tensile and bend- 
ing strength. Based on this information, it is possible to simulate the composite for these types of loading since the ma- 
terial continues to collapse after damage to the outer layer under bending and to the defective layer under tension. If 
there are no defects in the composite under tension, fracture mainly occurs in a dangerous section. 


Preparing for the experiment. For the experiment, samples were made from composite material (Fig. 1). 





b) 


Fig. 1. Test specimens: (a) bending; (b) tensile 
In the first case (see Fig. 1 a): 
- the type of laying is sandwich (symmetrical arrangement of layers in all samples), 
- the total number of layers per sample (first outer, then middle, then inner) are indicated, 
- 4 layers of basalt fabric (BF), 
- 2 layers of fiberglass (FG), 
- 7 layers of glass mat (GM). 
In the second case (see Fig. | 5): 
- 6 layers of basalt fabric (BF), 
- 5 layers of fiberglass (FG). 
In both cases, the binder is ED-20 resin, the hardener is PEPA. 
Thus, two types of material are obtained. Their differences allow you to check: 
- if using fabrics with mats is worthwhile, 


- what effects a combination of different types of materials will give. 
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Sampling. Samples were made through sequential placing of fabric layers (Fig. 2 a) on a base covered with a 


film, since there is no adhesion between the resin and the films. 





c) 


Fig. 2. Sampling: (a) material cutting; (b) removal of excess resin; (c) machine processing 


After placing all the layers, they were heated and smoothed at the same time with a spatula through the film 
(see Fig. 2 5); in this case, a smooth surface without excess resin was obtained. Sampling scheme: 

1) cutting samples with a margin for machining, 

2) preparation of the foundation, 

3) resin activation, 

4) placing layers, 

5) smoothing, 

6) fixation, 

7) curing, 

8) machining (see Fig. 2 c). 


Tensile and bending tests. Tensile and bending tests were carried out on a breaking machine (Fig. 3-5). 
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Fig. 3. Breaking machine 





b) 


Fig. 4. Testing: (a) bending; (5) breaking 
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b) 


Fig. 5. Results of destruction of samples: (a) tensile; (b) bending 


Research Results 

Test Conclusions. Based on the test results, the following conclusion can be drawn: when combining fabrics 
and mats, a characteristic feature of fracture is delamination at the border of various types of material (see Fig. 5). When 
in comes only to tissue layers, such an obvious feature is not observed. Tissues used as reinforcing agents are torn at the 
place of maximum stress and then are sharply destroyed. Initially, individual fibers are destroyed, which is accompanied 
by a characteristic sound. Then the fibers are pulled out of the matrix — and the matrix itself is destroyed. This is 
shown by decoloration — the place of destruction “turns white” and looks more matt. The test results are presented in 


Table 1. 

















Table | 
Characteristics of CM tested for bending and breaking 
6BT5ST/ED-20 4BT2ST7SM/ED-20 
Max. bending force, kg 500 965 
Max. breaking force, kg 3560 2600 
Max. bending stress, MPa 265.1 239.5 
Max. tensile strength, MPa 191.8 120 





Sample parameters and distance between sup- 


ports under bending, mm* 





b= 64.2;h=5.5;L=70 








b=61.7;h=8.2; L=70 








Section dimensions at break, mm 


b=35;h=5.2 


b=29.5;h=7.2 
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Density, kg/m° 





650 





700 








*Here, b is sample average width, / is sample average thickness, L is distance between supports. 





_ 
oO 


Composite material modeling and test validation. For modeling, the material 6BTS5ST was selected. Its 
characteristics are included in the Compass 3D library, and solid models are created. Using the APM FEM module for 
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Compass 3D, finite element calculations were performed. In this case, difficulties arose due to the selection of the opti- 
mal sizes of the finite elements. The following assumptions were used for calculations: the fixing and load points were 
performed by slight stretching of the “strips” of 0.1—0.2 mm wide. In these areas, the results of stresses and safety fac- 
tors should not be considered since they are incorrect. Such simplifications enable to apply loads and install fixation in 
any place on the models. The optimal mesh size is from 2 to 3 mm. 

The calculation results charts of stresses and safety factor of the ultimate strength are presented in Fig. 6. 


[Her 
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Fig. 6. Results charts: (a), (c) by stress; (b), (d) by safety factor 


Mechanics 


11 


http://vestnik.donstu.ru 


_ 
N 


Karnoub A., et al. Research and modeling of a multilayer composite material using basalt fabric 





The simulation results suggest that these samples are destroyed in places subjected to maximum stress. Due to 
the modeling error, the safety factor is other than unity, and the reduction in the size of the finite elements solves this 
problem. It is also worth saying that this method will be normally implemented when modeling structures operating in 
tension or compression, since these processes are accompanied mainly by normal stresses. However, during the tests, it 
has been found that at different elongation factors for the reinforcing materials, delamination occurs along the loading 
line at the layer boundary due to insufficient interlayer adhesion and high shear stresses. 

Discussion and Conclusions. The use of basalt fabrics as a reinforcing substance provides obtaining products 
with the properties of both glass and carbon plastics. Moreover, such a CM will be slightly more expensive than fiber- 
glass and much cheaper than carbon fiber. Good indicators of the specific strength of basalt, its absolute incombustibil- 
ity, high impact strength, and resistance to UV radiation should be mentioned. Basalt is a dielectric; therefore, it can be 
used in the manufacture of cases for radio equipment. 

The experiments have shown higher strength characteristics of those CM in which there are more layers of 
basalt fabric, and glass mat is not used. The use of glass mat as a core causes high interlayer shear stresses, which 
stimulates interlayer fracture both in tension and in bending. The maximum bending stresses for a sample with a large 
number of layers of basalt fabric amounted to 261.5 MPa (versus 239 MPa with a smaller number of layers of basalt 
fabric). In addition, there was a clear superiority in maximum tensile stresses — 190 MPa versus 120 MPa. 

The calculation results allow us to state that products made of composite materials (equated to isotropic mate- 
rials) can be modeled in CAD using the finite element method. It is very important to consider the type of loading on 
the product, since CM largely have anisotropic properties (the load is applied allowing for the direction of the fibers). 
So, in multilayer CM from structural fabrics, it is required to direct the loads along the fibers. In addition, it is necessary 
to consider the interlayer shear, different adhesion between the layers, etc. The main assumption of this method is the 
“constancy” of the thickness of the material, the number of layers and the order of their location. 

The data and method obtained need further research, which will allow us to create a range of composite mate- 


rials with known characteristics, to simulate them, and to perform the corresponding calculations. 
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Two ways of organizing scalar product in the boundary state method 


D. A. Ivanychev 
Lipetsk State Technical University (Lipetsk, Russian Federation) 





Introduction. The influence of two ways of organizing scalar product on the convergence rate of the solution in the en- 
ergy method of boundary states is considered. The method is based on the spaces of internal and boundary states which 
are conjugated through isomorphism. Both spaces are orthonormalized using one scalar product or another. The desired 
state is expanded in the Fourier series according to the elements of the orthonormalized basis; and the coefficients of 
this linear combination are determined. The two methods differ in the assignment of scalar products and the calculation 
of the Fourier coefficients. 

Materials and Methods. \n relation to the method of boundary states, a new theory of organizing a scalar product in the 
spaces of internal and boundary states is proposed. Computational algorithms are constructed for its practical imple- 
mentation. In the traditional (first) approach, the internal energy of elastic deformation is used as an orthogonalizer in 
the space of internal states. Here, the Fourier coefficients are the work of given forces on the basis vectors of displace- 
ment of the boundary points. In the studied (second) approach, scalar products are integrals of the cross products of the 
basis force vectors at the boundary. Accordingly, the Fourier coefficients are calculated as integrals of the product of 
the given forces at the body boundary by the basic force vectors. 

Results. A numerical study of the first primal axisymmetric problem of the elasticity theory for a transversely isotropic 
cylinder in the absence and presence of mass forces is conducted. In the absence of mass forces, an analysis of the elas- 
tic fields obtained for the same number of used basic elements has shown that the second method has the greatest accu- 
racy of the results. Under solving the problem with the presence of mass forces, the second method did not show effi- 
ciency in terms of the uniqueness of the solution; however, it is quite suitable for constructing a multitude of elastic 
fields used to solve more complex problems. 

Discussion and Conclusions. The results obtained can be used to solve boundary-value problems of mechanics of not 
only an anisotropic body, but also an isotropic one. When solving more complex problems, such as contact and mixed 


ones, the issue of the rate of convergence requires a separate study. 
Keywords: boundary state method, scalar product, internal energy, state spaces, the first main task, mass forces. 
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Introduction. The boundary-value problems of the elasticity theory in mechanics are sufficiently studied; 
therefore, in recent years, case studies have been carried out. For example, the first primal axisymmetric problem for a 
half-strip [1] is considered, the solution of which is constructed in the form of expansions in systems of Fadl-Papkovich 
functions and has an explicit form. A general method for solving the first primal problem of the elasticity theory for a 
rectilinear anisotropic body in the case of plane deformation is proposed [2]. Closed systems of boundary value prob- 
lems similar to Hilbert's problems are used, which makes it possible to achieve a more general method. Plane isotropic 


problems are solved using the finite element method based on the Castigliano variational principle [3]. This has provid- 
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ed for obtaining stress fields on grids of a sufficiently low dimension including those for incompressible materials. Con- 
tact problems on the implementation of elliptical stamps in a transversely isotropic elastic half-space are solved [4]. 

The method of boundary states under solving boundary value problems for anisotropic bodies also found its 
application. For example, the elastic equilibrium of a transversely isotropic cylinder under the action of axisymmetric 
surface forces is considered [5]. The problems of torsion of extended cylinders of a material with general anisotropy are 
studied [6]. A mathematical model is shown for obtaining explicit parametric solutions for isotropic and anisotropic 
bodies [7, 8], where the medium constants are included as parameters in the elastic fields. A technique has been devel- 
oped for solving problems of the elasticity theory using computer algebra [9]. In solving problems of a stress-strain un- 
bounded elastic medium containing spherical cavities or inclusions, the boundary-state method was used under different 
conditions [10]. 

In this paper, we study two approaches to the assignment of a scalar product in the “body” of the boundary 
state method. In this case, each state is tested using the example of solving the first primal problem of the elasticity the- 
ory. In each task, the same number of used elements is held and the error level is estimated. 

Materials and Methods. The boundary state method (BSM) [11] is an energy one; it uses the fundamental 
theory of series to solve the basic problems of mechanics. The concepts of internal and boundary states are used as sup- 
porting ones. 


The internal state € is due to a set of displacement vector u, strain tensor ¢ , and stress tensor T: 

b= (ue,T}. (1) 
The boundary state is determined by a set of displacement vector of the boundary points u" and forces p on the body 
boundary: 

y=’, ph. 
The totality of such states forms the basis of the spaces of internal == | Se ae eee Ae and boundary 
T= Gira tea Cents Freee states. Next, orthogonalization of state bases is carried out, where the following expression is 
used as an orthogonalizer in the basis of boundary states: 

(§.5,) se 

in the basis of boundary states — the expression: 


(1 1,) = [ure dS . (2) 


A single element €, € = corresponds to each element y, <€ J”, besides, this correspondence is one-to-one: €, > y,. 


This enables to reduce finding the internal state to the construction of a boundary state isomorphic to it. In the case of 


the first primal problem, the desired internal and boundary states are the Fourier series: 
G= 265.3 Y= Laie s (3) 
k=l k=l 
here, c, are the Fourier coefficients: 
c, =] pu,ds , (4) 
Ss 


where p is the vector of given surface forces; u, is the displacement vector in the k-th basic element of the basis of inter- 
nal states. In this case, basic sets are formed on the ground of a general or fundamental solution to the problem. 

The first general solutions to the Lame equation of the linear elasticity theory were built back in the 30s of the 
last century. The Lame equation is the Euler equation of motion (in this case, equilibrium): 


VT+f=0, 


where T is stress tensor; V is Hamilton operator acting as a divergence; f is mass forces. 
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In the Lame equation, the stress tensor T in accordance with Hooke's law is represented through the strain ten- 
sore. In turn, the strain tensor in accordance with the Cauchy relation is represented through the displacement vector u. 
In general solutions to the Lame equation, the displacement vector is determined: through the harmonic vector B and 
the harmonic scalar in the theory of isotropic elasticity, and through the stress function F — in the theory of anisotropic 
elasticity. 

Vector B (function F) can be represented as a series of basic vectors B, = B,(a') — coordinate functions a’. 
As a result, the following basic elements will be associated with each harmonic basis vector B, (function F, ): 

e displacement vector u, ; 

e strain tensor ¢, ; 

e stress tensor 7, ; 

e mass forces vector /, (from equilibrium equations); 

e surface forces vector (from Cauchy fundamental relationship): 

P, =n, » 

where n is the external unit normal to the surface of the body. 

According to the listed basic elements, the corresponding vectors or tensors are expanded in Fourier series with 


the same coefficients c, , which are determined from the conditions of orthogonality of the basic functions. For example, 


for the first primal problem in the absence of mass forces, when external forces p are given on the entire surface of the 
body S and orthogonalization of the basis vectors 


P, (J p,p,4S =, is Kronecker delta) is implemented, the coefficients c, are determined from the expression: 
Ss 


Cc, =Jpp,ds . (5) 
Ss 
This expression follows from the representation: 
P=XUCP, - 
k=l 


Thus, we study the solution generation method using the expression for scalar products in the basis of boundary states: 


(11) =| Pp dS (6) 


and the expression for Fourier coefficients (5). 


In the case of the second primal problem, there are the dependencies: 


(Yin) = [ujujas 5 
Ss 
c, = Juu;,ds , 
Ss 


where u is the given displacement vector of points of the body boundary; u; is the displacement vector in the k-th base 
element of the boundary states basis. 

Research Results. Parameters of the convergence rate of series and the result accuracy are considered using 
the example of solving the problem of the elastic equilibrium of a transversely isotropic cylinder of dark-gray siltstone 
[12] in a dimensionless form (Fig. 1). 

The boundary conditions are: 

p=9, S,| z=-2,0<r<l; 
p=9, S,| z=2,0<r<l; 
p, =4-z',p, =0, 8,| r=1,-2<2<2. 


ease 
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Fig. 1. Boundary conditions for the transtropic cylinder 


The methodology for solving the first primal problem in the traditional way using the scalar product (2) is de- 
scribed in detail in [6]. The basis of internal states (1) is constructed as follows: 

* using a general solution to the plane deformation problem [13], basis sets of plane auxiliary states are con- 
structed; 

* basis sets of spatial axisymmetric states are determined from the transition formulas; 

* orthonormalization of the internal states basis is carried out according to the Gram-Schmidt matrix algorithm 
using the scalar product (2); 

* the orthonormal basis of boundary states is reduced from the orthonormal basis of internal states; 

¢ Fourier coefficients (4) are calculated, and series (3) are constructed in an expanded form (index k is placed 


on top): 
< k < k < k < k 
U, = Leh, > DP; = LP 5 0, = VC, 0; 5 €y = VCE; - 


We omit the information on the fields of the stress-strain state characteristics obtained with one and the other 
method of assigning the scalar product, and we give only the main results. In addition, we call the traditional approach 
used in [11] as the first method for solving the problem, and the approach using the scalar product (6) and Fourier coef- 
ficients (5) as the second method. The solution accuracy while keeping the same number of basic elements is higher in 
the second method. In Fig. 2, a comparison of the obtained boundary conditions (BC) with the specified ones when us- 
ing the 8 values of the Fourier coefficient is given for each method. The efforts are shown to scale, for example, the true 


value p, in the first graph of Fig. 2 is equal to the value on the graph multiplied by the coefficient « . 
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Fig. 2. Verification of BC for cylinder with 8 basis elements 
This trend is preserved with an increase in the number of basis elements used. For the 61st element, verifica- 
tion of the boundary conditions is presented in Fig. 3 (comparison of the force p, on the border section S, is given). If 
we evaluate an error as the maximum deviation of the obtained value from the given value, then in the second method 


the error is less. 
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Fig. 3. Verification of BC for cylinder with 61-st basis element 


Next, we examine the solution accuracy for a transversely isotropic axisymmetric body of noncanonical shape 


(Fig. 4). Border conditions are: 


pa0, 8 95;5 
p.=0,7 =025,5,| 2=1,04r21; 


Pp, =9,p, =-1, S,| z=-1,0<r<0,5. 





Fig. 4. Boundary conditions for axisymmetric body 


15 elements of the basis are kept. Fig. 5 presents a comparison of the boundary conditions for each method 


(not all boundary sections and components of the force vector are shown). 
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Fig. 5. Verification of BC for axisymmetric body 


As can be seen from the graphs, for a body of a more complex shape, a difference in the convergence rate is 
observed in favor of the second method. Consider tasks involving mass forces. The sequence of the solution generation 
is as follows: 


1. The dependence of the displacement vector of the planar auxiliary state on the coordinates y“z? is specified 


and, on its basis, the displacement vector of the spatial axisymmetric state is determined. 
2. For sucha vector, the following parameters are determined: 
* strain tensor according to the Cauchy relation; 
¢ stress tensor from Hooke's law; 
* efforts on the body surface from the fundamental Cauchy relation; 
« mass forces from the equilibrium equation. 
3. An exact particular solution to the problem corresponding to the displacement function specified at each point 
of the body is constructed. 


4. Sortinga+B<n (n=1, 2,3 ...), a set of exact particular solutions of the problem of the linear theory of elas- 


ticity for the parameters is constructed: 
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* displacement vectors; 

* strain tensors; 

¢ stress tensors; 

* vectors of surface forces; 


* vectors of mass forces. 


5. The bases of internal Be cas eet and boundary Pi ea ot states in which the 


equalities are respected, are generated: 
E, = {uy & Ty}, 


Yer {Uj Pio X ih . 
6. We leave only linearly independent among these solutions and carry out their orthogonalization in accordance 


with scalar products in the bases of internal and boundary states: 


(656)) = J Eo dV , 
(¥,) = J pe uPads + | XPuPav 
S V 


(indices i andj, responsible for the numbers of elements, are placed on top and are enclosed in brackets). 
7. As aresult, we obtain a basis through which the corresponding vectors or tensors are expanded to series (3) 
with equal coefficients: 
Cc, =[/Xu,dV +) puzds , 
V S 


where X is the vector of given mass forces. 

We study the possibility of constructing an elastic field in the presence of mass forces using relations (6) and 
(5). Consider the first primal problem with unbalanced forces for a transversely isotropic cylinder (Fig. 6). Border con- 
ditions are: 


p=9,S,US;; 


p= p.=7 S| = 2, 05721, 





Fig. 6. Boundary conditions for cylinder 


The expressions for the orthonormal basis set of the displacement vector components u = {u,w} are given in 


Table 1. 
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Table 1 
Orthonormal basis set of displacement vector components 
u WwW 
oh 0 0.1213z 
é 0.0788r ~0.1126z 
53 0 0.033127 
E, 0.0592rz —0.0268z? 
Ss —0.02rz 0.0975r? + 0.0092” 
Ss 0.15357 ~0.9163z + 0.0691z° 
G7 —0.2896r + 0.072412" 1.2067z—0.1005z° 
Es 1.86297 —0.4824rz* —8.2447z + 0.4824r*z +0.6701z° 
E, 0.09087rz —0.0454r? — 0.434727 —0.051 1z* 
E10 —0.1757rz + 0.0439rz° —0.0195r? + 0.3661z* —0.0457z* 
Gi 0.288 1rz —0.0804rz* —0.144r? — 0.72092? + 0.1206r7z* +0.0837z* 
ye 


To obtain a rigorous solution, it took 11 Fourier coefficients, non-zero values: c, =0.2178; c, =—0.1226; 
c,; = 0.2377 ; c, =—0.0732 ; c, =0.0247; c, =0.1443; c,, = 0.1443. We give the expressions for displacements and 
mass forces (4 decimal places are kept): 
u =0.2592r + 0.0367 rz —0.0696rz* —0.011672’; 
w=-0.018377 —1.1497z +. 0.0696 r7z —0.094z7 +; 
+0.01747r°z* +.0.0967z° + 0.01212"; 


R=-0.2814r—0.1407 rz; Z =1.2012—0.257r° —3.6038z—0.9z’. 


If we construct a basis of internal states through planar auxiliary states formed using monomial z“y’, then the 


orthonormal basis and the Fourier coefficients change. In this case, the decision will also be rigorous and will be ac- 
cepted in: 
u = 0.0363r —0.0913 rz + 0.0309 r°z —0.0287 12’ ; 
w=0.0456r° —0.0077 r* —0.4148z + 0.0287 rz +0.0257z* + 0.03992’ ; 
R=-2.0145r-1.61787z; Z =—0.25r° -1.4893z. 

Similarly, other particular solutions to the problem can be obtained using different types of polynomials, for 
example, z“y’ +.z%, etc. to form the basis. 

The approach under study provides for lots of solutions of one boundary value problem of the elasticity theory 
in the presence of mass forces. Displacements and mass forces whose combinations give the distribution of stresses 
satisfying the given forces at the boundary are also subject. The first method is devoid of this feature because mass 
forces in it are part of the given conditions, and the task is only to find the field of displacements. 

Discussion and Conclusions. The second method to solve the problem has the best convergence. In addition, 
unlike the first method, when calculating scalar products in the process of orthogonalization and in determining the 
Fourier coefficients, the second method does not use deformations and displacements. Here, a basic set of stresses is 
generated, and its trace at the boundary is a basis set of surface forces with which orthogonalization and series construc- 
tion are performed. This means that when calculating scalar products, there is no error associated with the components 
responsible for the rigid displacement that may occur during the formation of the basis [11]. 

If problems with mass forces are considered, then the second method may be useful in the formation of many 
particular solutions whose stresses satisfy certain conditions at the boundary. These solutions can be used as the basis 
ones for a more complex problem, and also be useful under determining the elastic fields realized from fictitious loads 
resulting from the application of the Poincare method [7, 8]. 

The accuracy of solving problems of the elasticity theory by the boundary-state method is analyzed using dif- 
ferent approaches to the construction of scalar products. The solution to the problem of the linear elasticity theory using 


Mechanics 


23 


http://vestnik.donstu.ru 


Ss) 
ws 


Vestnik of Don State Technical University. 2020. Vol. 20, no. 1, pp. 15-24. ISSN 1992-5980 eISSN 1992-6006 





the representation of a general solution to the Lame equation in the form of a Fourier series on basis functions and ex- 
pression (6) as the orthogonalizer of these functions had the best convergence. 
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Introduction. The history of solving the problem of welding structures made of stainless and heat-resistant metals and 
alloys goes back several decades. Researchers were particularly interested in working with deformed aluminum alloys 2— 
6 mm thick. As a rule, such thin-walled structures are welded in an argon shielding gas at relatively small currents; there- 
fore, metal transfer is large-droplet (the weld is shaped in the form of separate large droplets with a narrow penetration of 
the welding components). At the same time, the weld is very convex, which does not meet the operational requirements 
of the structures. Thus, it was important to solve the following problems: to obtain a controlled fine-drop transfer of elec- 
trode metal at currents corresponding to a large-drop transfer; to determine a condition for the controlled transfer; to de- 
velop a power supply system for the welding arc. 

Materials and Methods. The behavior and parameters of the arc were recorded through the high-speed film and video 
shooting with synchronous oscillography of the electrical process parameters — current and voltage. They were recorded 
by light-beam oscilloscopes and two-screen oscilloscopes. The data were processed using a computer complex and Dia- 
dem 10.1 software. 

Results. The basic condition for the controlled metal transfer is determined through applying current pulses to the welding 
arc from special pulse sources with and without energy storage devices. Transients in the electrical circuits of the main 
welding source during the current pulse action and pause are considered. The factors providing stability of rigid and flex- 
ible pulsed-arc welding (PAW) are indicated. 

Discussion and Conclusions. The results of studying the possibility to control the welding arc processing behavior and the 
proposed methods for calculating the parameters of the PAW mode became the basis for the development of technology 
and equipment for the mechanized GMAWP of aluminum alloy assemblies. They are introduced at the enterprises of the 
aviation industry, shipbuilding. Solutions for stainless and heat-resistant steels and alloys are used at the motor industry 
enterprises. As a result of studies on the mechanized CO, activated electrode wire welding, a mechanized PAW technolo- 
gy was developed for units of stainless steel electric furnace bodies, structures of road-building, agricultural machinery 


and ships. 
Keywords: welding in shielding gases, consumable electrode, current pulse parameters, controlled metal transfer. 
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Introduction. Since early sixties of the XX century, gas-shielded consumable-electrode welding of metal struc- 
tures has been studied [1]. One of the basic problems of this process is the manufacture of welded structures from stain- 
less and heat-resistant metals and alloys. Of particular difficulty was the creation of deformed aluminum alloys 2-6 mm 
thick. These materials are used in aviation, motor and shipbuilding, in light, food, and chemical industries. Welding of 
thin-walled structures of these metals was performed, as a rule, by a non-consumable electrode in the argon shielding gas 


medium with relatively low process productivity. Mechanized consumable-electrode welding at relatively small currents 
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was also used, which means with large-drop metal transfer. In this case, the weld at small currents is formed in the form 
of separate fused large droplets with a narrow, small penetration of the welded elements, which does not always meet the 
technical requirements for the structure design. 

It is known that there is a critical current under welding in argon. In a narrow range of its changes, the character 
of metal transfer sharply changes. At currents below the critical level, the metal transfer is globular, and above the critical 
one, it is spray transfer. In case of the globular metal transfer, a weld of improper shape is obtained — with a narrow, shal- 
low penetration and a large convexity. If the current is higher than critical, during spray transfer, the weld has finger- 
shaped penetration, which reduces the durability of the weld joint. 

During the mechanized welding of a melting electrode wire, the current value periodically changed for a short 
time. The purpose of this operation is to obtain controlled fine-droplet transfer of electrode metal at currents corre- 
sponding to globular metal transfer. In this case, the geometric shape of the weld should meet the technical require- 


ments for welded structures. 


Materials and Methods 

The arc power system for pulsed arc welding (PAW). To achieve this goal, an arc power system from two 
electric energy sources has been developed. The system includes: 

- the normal power supply is the PSG-500 welding transducer with a flat-dipping volt-ampere characteristic 
(VAC) since at that time there were no welding rectifiers for the mechanized welding in shielding gases; 

- a special switching power supply (SPS) is a pulse energy storage generator. 

The normal and pulsed sources were connected in parallel to the arc gap. 

In the first arc SPS developed, 150 A current ignitrons were used as power controlled rectifiers. Then, when 
the VK-200 and VT-150 power semiconductor rectifiers became available, they began to manufacture and introduce 
switching power supplies for semiconductor valves at enterprises. Switched energy storage power supplies are protected 
by copyright". 

When the normal and SPS of the arc are connected in parallel [2], an increased voltage is applied to the arc un- 
der the current pulse pile-up. It is directed opposite to the normal source and “locks” it due to the presence of an uncon- 
trolled rectifier V2 in its circuit (Fig. 1). A decrease in current in the normal source circuit under the action of a current 


pulse is meant by “locking”. 





Fig. 1. Equivalent circuit of arc supply from normal and switched power sources: Egan, Ev, £, are voltage of welding, switched 
power sources and arc; igun, tuun» i, are current of welding, pulsed sources and arc; R.,, R1, R, are active resistance of the welding 
circuit, pulse source and arc; £2, L1 are inductance of the circuit of the welding and pulse sources; V2, V1 are rectifiers in the circuit 
of the welding and pulse sources 


However, despite the presence of the V2 rectifier, the “locking” of the normal source does not occur instantly. 
Due to the inductive resistance L2, a transient process occurs in the circuit (during the pulse pile-up and after its action). 
Consider the transition process of current variation in the normal source circuit when a current pulse is applied. 


To do this, we solve a linear differential equation compiled for the normal power supply circuit (see Fig. 1): 


' Budnik NM, Dyurgerov NG, Sagirov KbN, et al. Ustroistvo dlya impul'sno-dugovoi svarki: A.S. 226752 SSSR [Device for pulsed-arc welding: au- 
thor’s cert. no. 226752] USSR, 1968. 

° Budnik NM, Sagirov KhN, Dyurgerov NG, et al. Ustroistvo dlya impul'sno-dugovoi svarki: A.S. 299111 SSSR [Device for pulsed-are welding: au- 
thor’s cert. no. 226752] USSR, 1971. 
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BE, = by 22+ Role + Radu: (1) 


Here, Ey = Eoyy — Ey, Ro = Ryt+Rey; Ig is current of the normal power source; iyy, is free current in the pulse 


(for example, for sinusoidal pulses with a damping amplitude): 
inun = ar exp (- Sa t) sinWot, 


_ - }4 Ray 
where Ry, = Ry, + Ry, @o= ae et 





After the pulse end, the current in the normal source circuit increases exponentially, it is determined by the 


time constant of the normal power source circuit Ta, = Rey/L2. 


Is (tt, “ty = [Ino — Is (t)e=t, | (1 — exp — +) + Ig ()t=t5- (2) 

Based on the equations (1 and 2), the law of the current variation in the circuit of the normal source was ob- 
tained during the current pulse pile-up and after its end [3, 4]. 

Research methodology. The arc behavior and its parameters under welding at various spatial positions of non- 
ferrous and ferrous metal samples were recorded by high-speed film shooting up to 5000 frames/s and video shooting 
2000 frames/s. At the same time, synchronous oscillographic testing of the process electrical parameters was performed: 
the arc current and voltage were recorded by the light-beam and electronic oscilloscopes. The data obtained by means of 
a two-screen electronic oscilloscope was processed using a computer complex and Diadem 10.1 software. 

We studied the nature of metal transfer, the weld formation, and determined the PAW physical and process pa- 
rameters in the protective gases of these metals. 

Research Results 

Welding in argon protective environment. An analysis of the obtained patterns provides evaluating of the ef- 
fect of the PAW mode parameters on the joint operation of the arc feeding system, in which the welding process be- 
comes unstable. Oscillograms of the PAW process at various inductances of the welding power supply circuits are 


shown in Fig. 2. 


— — ~aP gn > "ar 4 
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Fig. 2. Oscillogram of current (/,), voltage (U,) of the are and current in the circuit of the welding power source (/,,) (t, = 4.6°10°°s, 
I, = 690 A, f= 100 cps): continuous arcing J.) = 70 A, L = (1.0 = 1.1)-10° H (a); intermittent are burning Tgp = 150 A, 
L=0.3-:10-*H (b) 

The natural (non-pulsed) and pulsed welding processes are considered. It is established that for identical values 
of the effective current, the velocity of droplet movement in the arc gap is always higher under PAW. Moreover, a larg- 
er current in the pulse at the moment of its detachment from the electrode always corresponds to a larger drop velocity. 
Separation of the droplet from the electrode can occur in various phases. It depends on the welding mode and current 
pulse parameters. 

When the droplet is detached at the end of the pulses, the current is close to the base current of the welding 
process, and the droplet flying speed reaches 1.2—2.0 m/s. In this case, a controlled directional metal transfer is provided 


under welding in all spatial positions with minimal metal loss through splatter (Fig. 3 a). 
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Fig. 3. PAW records of aluminum alloy AMgé6 in the overhead position under various modes (shooting speed is 4000 frames/s, 
welding wire is AMg61 1.6 mm in diameter): effective current is I,4= 100 A, t= 2.5:10°% s, I, =310 A, f = 50 eps (a); I,m = 180 A, 
ty= 4.510, I, =500 A, f= 50 eps (d) 

It is established that the minimum amplitude of the current pulse [, min, which provides a transition from un- 
controlled transfer to controlled transfer, depends on: 

— surface tension of the electrode metal o,, 

— electrode diameter d,, 

— welding mode parameters, /,», v 


[2 


Se = Auge fe (3) 

In the case of droplet detachment from the ce at an amplitude value of the pulse current, its flying speed 
in the arc gap can exceed 8 m/s. This causes the appearance of: 

— near-weld splashes (due to spattering of molten metal from the weld pool), 

— a weld-affected zone of fine droplet spraying formed when a neck of a molted metal breaks between an 
electrode and a droplet (see Fig.3 5). 

To reduce the metal loss due to spatter, the amplitude of the current pulse should not exceed 1.1 J, min. 

The increase in [,, is limited by the condition of continuous arc burning: if the effective current is constant, 
then with an increase in the amplitude, frequency, and pulsewidth, the base current decreases. 

When the base current decreases below the minimum steady-state current of the arc, its breaks occur (the weld- 
ing current is interrupted), the stability of the welding process is violated (see Fig. 2 5). This condition was used to de- 


termine the maximum values of the pulse parameters. For example, the maximum amplitude of a pulse from a source 


with energy storage for unipolar current pulses (damped sinusoid): 


Ic (= 
Linge < ae [(nt2 + ty W/ fen (1 — exp ( zit)). (4) 


Fig. 4 shows the dependences of the minimum required amplitude of the pulse current (curves 1, 1') and the 





maximum allowable amplitude (curves 2, 2') on their repetition rate for different pulsewidth. Curves 1, 1' (with increas- 
ing current) are the boundary of the transition from uncontrolled transfer of metal droplets to the controlled transfer. 
Curves 2, 2'are the boundary of the transition from controlled droplet transfer to arc burning cessation. These curves 


limit the range of the PAW modes by a consumable electrode. 
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Fig. 4. Areas of permissible amplitude values of current pulses: argon, 
alloy AMg6, Ing = 140 A; 1, 2 —t,= 1.4 ms; 1,2 —t,=2.2 ms 

From the equations (3) and (4), it follows that with an increase in the pulse repetition rate, the minimum re- 
quired amplitude of the current increases, and the maximum allowable amplitude decreases. The area of controlled drop 
transfer narrows. 

Under the limiting current pulse rate (f,p), the frequency is adopted at which the minimum required amplitude 
of the current pulse is the maximum allowed. This frequency was determined through solving equations describing 
In min (3) and Inmax (4). 

Welding in active shielding gases. Welding in active shielding gases (CO:, N, and mixtures of 
Ar < 80 % + CO, > 20%) and He on the current of direct and reverse polarity is considered. It was established that in 
this case, the current pulse pile-up with the parameters used for welding in argon does not provide the controlled metal 
transfer [4, 5]. The current pulse pile-up under the activated electrode wire welding (alkali salts) by the current of direct 
polarity provides stabilization of the arc — it takes the shape corresponding to the current higher than critical. However, 
in the pulse spacing, when the arc current is less than critical, it wanders intensively. If a drop passes from the electrode 
at the end of the pulse, when the current decreases to below a critical value, then the detachment and transfer of the drop 
takes place in an erratic arc. In the pulse spacing, the electrode wire melts, and at its end, there is always a certain 
amount of molten metal. Due to the wandering of the arc, the drop moves relative to the axis of the electrode wire. 
When a pulse pile-up, the drop detaches and may not fall into the weld pool. The pulse pile-up with duration of 1.2— 
1.8 ms causes the detachment of droplets from the electrode wire with each current pulse; however, it does not provide 
directional transfer of metal into the weld pool and is characterized by increased spatter. 

It was found that non-coaxial droplets, under current pulse pile-up, deviate to a greater extent from the elec- 
trode wire axis than coaxial ones. The largest number of misaligned droplets detaches from the electrode at an angle of 
20-30° and does not fall into the weld pool. Coaxial drops detach at an angle not exceeding 10°. 

When the pulsewidth is 4.0—-5.0 ms, the base current of the welding process decreases due to the integral self- 
regulation of electrode melting [6, 7]. At the end of the electrode wire, despite the wandering of the arc during cessa- 
tion, a small volume of molten metal is formed. As a rule, under the current pulse action, several coaxial droplets are 
formed, which are transferred to the weld pool, and this causes small spatter. 

To reduce the metal loss due to spatter under the activated electrode welding in active shielding gases and their 
mixtures, it is required to combine the processes of melting and droplet transfer into the weld pool. The electrode wire 
should melt at a current above critical, when the arc is spatially stable and does not cause deviation of the flight trajecto- 
ry of droplets from the electrode axis in the interelectrode gap [8, 9]. This is achieved as follows: current pulses of the 
same polarity are superimposed on a continuously burning arc with the minimum possible amplitude J, = (1.5—2.0) ep 
and the maximum possible duration ¢,, = 4.0—-10 ms at a pulse repetition rate of 100— 50 cps. The base current is selected 
to be minimal (so that during cessation, the electrode wire does not practically melt). In this case, during the period of 
the pulse action, short-term spray transfer (ST), called intermittent spray transfer (IST), occurs. When the pulse ampli- 


tude decreases below 1.5 J, 


xp» metal spray transfer is impossible. 


A nonstorage switched power supply has been developed and manufactured for the practical implementation of 


the described process. 
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For the first time, the features of the weld formation under PAW were determined. The nature of the molten 
metal displacement in the weld pool was studied with and without current pulse pile-up on the arc. This provided de- 
termining factors that improve the weld formation and the retention of the metal of the weld pool under PAW in various 
spatial positions. 

Using high-speed filming (up to 5000 frames per second), it was found that under PAW, molten metal is dis- 
placed from the head end of the weld pool under the action of arc pressure and pressure created by drops of a trans- 
ferred electrode metal [3]. It turned out that the speeds of molten metal movement caused by the arc pressure under the 
action of the current pulse and the kinetic energy of the droplets were commensurate. Therefore, in case of PAW, in one 
cycle, the molten metal is twice displaced from the head end of the pool to its tail end. 

The move frequency of molten metal in the pool increases with an increase in the pulse frequency, ceteris pari- 
bus (constant welding current and arc voltage). This cuts cycle time and reduces the amount of transported metal, which 
crystallizes in a thinner layer, which helps to improve the primary crystallization of the weld metal. 

PAW when feeding an arc from two sources (see Fig. 1), as a rule, is performed according to a rigid program 
with a fixed pulse repetition rate. The shape of the current pulses is similar to the curve of the applied increased voltage 
(see Fig. 2). 

Since 1989, research has been carried out to create welding rectifiers with inverter converters for various with 
consumable and non-consumable electrode welding methods!”. 

PAW by pulses of a rectangular shape. Welding rectifiers with an inverter converter are created on the basis 
of high-frequency inverters with a current conversion frequency of 16 kHz. In this case, power sources of the welding 
arc are low-inertia — high-speed ones. This provides a curve of the welding current of various shapes with a change in 
magnitude over a wide range. The rectifiers with combined volt-ampere static characteristics with three prominent sec- 
tions are most used when performing welding operations. The steep-dipping section at currents less than 30A provides 
reliable arc striking and the establishment of the welding process. The flat-dipping section and the bayonet section in 
the range of operating modes provide stabilization of the welding current in a wide range of its regulation, in a wide 
range of changes in the voltage drop across the arc. 

The digital control system for the arc power source provides a fast shape change of the welding current curve 
and maintaining the bayonet section of the volt-ampere characteristics in the range of the current control limit. The 
source is equipped with a control system for the consumable and non-consumable electrode welding. 

PAW by rectangular pulses from a single arc power source can be obtained for the mechanized consumable- 
electrode wire welding with a constant feed rate. For this, it is required to use a discrete two-level switching of the volt- 
ampere characteristic of the power supply (PS VAC) with a bayonet section in the range of operating currents of the 
welding rectifier. 

Let us consider how the welding electrical parameters change at a constant feed rate of the electrode wire and 


discrete switching of the VAC of a welding rectifier with an inverter converter (Fig.5). 


'Lenivkin VA, Petrov PI, Oleinikov AG. Sposob dugovoi svarki plavleniem: A. S. 1776517 SSSR: [Fusion arc welding method: author’s cert. no. 
1776517] USSR, 1991. 

? Oleinikov AG, Lenivkin VA, Petrov PI, Zhmailov BB. Istochnik pitaniya dlya dugovoi svarki: RF patent 2063850 [Arc welding power source] RF 
Patent no. 2063850, 1994. 


Lenivkin V. A., et al. Study on pulsed-arc welding issues at the Machines and Welding Production Automation Department, RIAE — DSTU 





Is In I 
Fig. 5. Variation of mode energy parameters under PAW according to flexible program: | and 1' are PS VAC providing base current 
J; and pulse current J/,, respectively; 2 is the given curve of self-regulation characteristics of electrode wire melting; 2' and 2" are 
dummy curves of the self-regulating melting behavior of the electrode wire; 3 and 3, are static are VAC 
before and after the current pulse action 


The electrical parameters (voltage U, and arc current /;) of a stable mechanized consumable-electrode wire 
welding are determined by the point a (here, PS VAC I and the characteristic line of self-regulating melting of the elec- 
trode wire 2 meet at a constant feed rate). 

It is known that at any point of stable welding operation, the electrode wire feed speed (v,,) is equal to its melt- 
ing rate (v,), and the arc gap length (€,) depends on the voltage drop across the arc and remains unchanged. Therefore, the 
static volt-ampere characteristic of the arc (A VAC) 3 with its constant length ¢, passes through quiescent point a (Fig. 5). 

The current value at point a (see Fig. 5) is selected less than the critical current value, that is, the process of arc 
welding with a globular droplet metal transfer for a specific diameter and type of electrode wire corresponds to the base 
current of the welding process J;, which is determined approximately through the dependence: 

15 = Vy / Key (5) 


1 Onn 


Here, v, is the electrode feed rate; k,, is the coefficient of arc current self-regulation; k.. = Fay ae 
" 3 


(Onn is the 


electrode wire fusion coefficient, g/A h); y is the electrode wire density, g/cm’. 

The voltage drop across the arc (U,) at point a is selected so that PS VAC I can provide stable arc burning 
without short circuits of the arc gap. 

Under a discrete switching of PS VAC from position | to position 1' and a constant arc length, the welding 
mode parameters of point a move to the position of point a,. (In Fig. 5, the arrows indicate the direction of change in 
the magnitude of J, and U,). The arc current value at the point a, corresponds to the given amplitude value of the pulse 
current (/,,), which is selected from the condition for providing the ST under its operation: 


Iq 2(.5£2.0)1,,, 


where J,, is the critical current, which is determined by the welding current magnitude, the diameter and grade 
of the electrode wire, the length of the arc gap. 

A stable welding process with continuous arc burning at point a, and an unchanged arc length f, and /,, can be 
established at the intersection of PS VAC 1' with the self-regulating melting line of the electrode wire 2' if the voltage 
drop is slightly greater than at point a (v, = v;). This is due to the fact that under real conditions of welding, when meas- 
uring the voltage drop across the arc, the voltage drop at the electrode wire extension is simultaneously taken into ac- 
count. Therefore, an increase in the voltage drop at point a, is caused by an increase in the voltage drop at the electrode 
extension and the arc resistance at a higher current /,. 

When PS VAC are switched from position | to 1', the welding mode parameters (arc current and voltage) are 
determined through the intersection of PS VAC 1' and A VAC 3. In this case, the welding process is unstable since the 
melting rate of the electrode wire at point a, remains unchanged and is determined by the self-regulating melting curve 


of electrode 2, and v, at J, will be more than v, by the value Av,: 
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Av, =k.,, 0, —1,)- (6) 

At that, the arc length will start to increase. The static characteristic of arc 3 will start to shift equidistantly up 
from point a, in the direction of point 6). 

It is established that with an increase in the arc length at a constant value of the welding current and a constant 
feed rate of the electrode wire, its stick out decreases, and the thermal power spent on melting the electrode decreases. 
The frequency of metal spray transfer decreases and the size of the transferred droplets increases. Spray transfer trans- 
forms monotonously into globular one. 

When combining point a, with the intersection point of PS VAC 1''and the self-regulating melting line of elec- 
trode wire 2, a new point of stable operation with metal transfer is established, as at point a, at which v, = v3, but with a 
different arc length. 

For PAW with ST, it is required to switch the PS VAC from position 1' to position | with a shorter arc length 
é,3 than at point 6), at which the ST is not violated yet. 

During the arc burning on a current J, A VAC moves equidistantly from point a; to point 6). In this case, the 
arc length increases by Af,,, and the arc static characteristic will occupy position 3, (Fig. 5). The voltage drop across the 
arc will increase by 

AU, 7 k Ata Koen = ig); (7) 
where k, is the arc core voltage gradient, V/mm; k,, is the coefficient characterizing the slope (A VAC). 

Under a discrete switching of PS VAC from position 1' to position 1, the mode parameters and the welding 
process go from point 6; corresponding to the pulse current /,, to point 6 corresponding to the base current J;. Welding 
at point 6 would be stable if PS VAC intersected by the line of the self-regulation characteristic of the electrode wire 
melting 2, which corresponds to lower v, than the specified curve speed of the self-regulation characteristic of electrode 
wire melting 2. 

Therefore, the melting rate of the electrode wire v, at point 6 slows down sharply at almost constant base cur- 
rent J;. The arc gap starts to contract since v, > v3. The arc static characteristic is shifted equidistantly from point 6 to 
quiescent point a. During the arc burning at the base current, the arc length decreases by Af,p. 

When combining the static A VAC 3, with characteristic 3 Al.) = | Alo ls the voltage on the arc becomes 
equal to the specified lower level of operation of the two-level sensor. At this moment, PS VAC switches from position 
1 to position 1', and the process is periodically repeated. 

The pause time t,, is determined by the burning time of the arc between points 6 and a at current J;, and the 
pulse length t,, — between the points a, and 6, at current J;. 

PS VAC from position | to position 1' and back, switches automatically depending on the specified voltage 
drop of the arc burning at the base and pulse current. The voltage drop is set for the base current a and pulse current 6. 
The power supply control system receives data on the base and pulse voltages using a two-level voltage sensor. In this 
case, the control system of the welding power source operates according to a flexible program. The switching frequency 
«pulse — pause» depends on the set difference in the voltage drop between the values of the base and pulse currents. 

In some cases, for example, when welding in the vertical and overhead positions, a higher flight speed of drop- 
lets is required than under the natural metal transfer. Therefore, the pulsed process should also be carried out in case 
when the welding current corresponds to the fine-droplet metal transfer and the following condition should be satisfied: 

ee lh 21s (8) 
where k,,, is the coefficient of droplet transfer irregularity, f, is the droplet transition frequency under the impact of 
current pulses, fis the droplet transition frequency of a natural process. 

The optimal parameters of the pulse process, depending on the spatial position of the weld, should meet the 
conditions described by the equations (3) and (8). 

The minimum pulse length providing intermittent spray transfer is (4.5-5.0) 10° s at J, = (1.5-.0) Ip. In this 


case, the pulse repetition rate is maximum (110-100 cps). The maximum pulse length is determined from the condition 
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for providing ST at a constant value of J, and the maximum allowable ¢,. The pulse repetition rate is selected from the 
condition of uniformity in the formation of the weld width along its length and should not be less than 10 cps. 


Fig. 6 shows fragments of the PAW current and voltage oscillographs with intermittent spray transfer’. 
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Fig. 6. Section fragment of synchronous oscillograms of PAW by rectangular current pulses with IST: 
(a) voltage, (b) arc current 
On the oscillograms, the moment of transition of arc burning from the base current of 160A (point 1) to the 


pulse current (/, = 375A, point 2) is recorded. The current of 160A is less than the critical current for an electrode wire 
of 1.2 mm (J, is 190+200A). In this case, the current in the pulse is (1.97+1.87) J. The voltage drop across the arc at 
the time of switching from J, to J, is 22.8V. When the arc burns during the current pulse action, the voltage across it 
increases to 25.8V. 

Therefore, during the action of the current pulse t,, with a duration of 31-10~ s, the are length increases by 
At, <2 mm, and during the arc burning by J; in a pause t, with a duration of 29-10™ s, the arc length decreases by the 
same amount. The pulse repetition rate is 16,6 cps. 

PS VAC switches automatically from position 1, corresponding to welding at a pause current (/g), to position 
1, (corresponds to the pulse current) and back. The switching depends on the set voltage drop 

AU, =k, Aly 
where k, is the voltage gradient of the arc core under welding in argon at a current above critical 1.4 V/mm) between 
points a and 6, (corresponds to the arc length at these points, see Fig. 5). 

Using a two-level voltage sensor, the power supply control system receives information on the variation in the 
voltage between the base and pulse currents. 

Discussion and Conclusions. Two main varieties of PAW are known. The first one is a “pulse-drop” opera- 
tion: at t,, equal to (1.5—2.5):10~ s, each pulse at the end of its action transfers one drop of electrode metal from the 
electrode wire to the weld pool. In this case, the processes of melting and transfer of the electrode metal are largely sep- 
arated in time. 

In the second type of PAW, the pulse acts much longer (t, > 4.0-10~* s and more), and the intermittent spray 
metal transfer occurs. In this case, the processes of melting and metal transfer are combined. 

The control system of the pulse-arc welding process provides a smooth adjustment of the repetition rate of cur- 
rent pulses, which depends on the specified increment of the arc gap length, i.e., on the operating voltage of the two- 
level sensor and the amplitude of the pulse current. The repetition rate of current pulses, the duration of pulses and 
pauses change smoothly and are determined by incrementing the arc gap length under the action of the pulse current 


depending on the values of the specified current in the pause and pulse. 


' Lenivkin VA, Kiseley DM, Dyurgerov NG. Sposob impul'sno-dugovoi svarki: RF patent 2570145 [Pulse arc welding method]. RF Patent no. 
2570145, 2015. 
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The flexibility of the IST program is achieved through the use of a two-level voltage sensor. The lower level of 
the sensor response voltage is determined by the parameters of the welding mode, which provide a stable process in a 
pause at a current below the critical level. The upper level depends on the elongation of the arc by 0.5—2.0 mm under 
the arc burning at a current that provides ST. 

The developed PAW system with rectangular pulses through switching the integrated VAC-rectifier with an 
inverter converter is used for the mechanized welding in argon and a gas mixture (argon> 80% and CO, <20%) accord- 
ing to the rigid “pulse-drop” program and the flexible “intermittent spray transfer” program. 

The study results of the capability of controlling the processing behavior of the welding arc and the proposed 
methods for calculating the PAW mode parameters laid the groundwork for the development of technology and equip- 
ment for the mechanized consumable-electrode PAW of aluminum alloys AMts, AMtsM, AD1, AMg6, AMg61, 01915. 
These developments were introduced at the enterprises of the aviation industry, shipbuilding. Solutions for stainless and 
heat-resistant steels and alloys are used at the motor industry enterprises. 

As a result of the studies on the features of the mechanized activated electrode wire welding in carbon dioxide, 
the mechanized activated electrode wire Sv-O8G2S PAW technology was developed for the electric furnace body units 
made of stainless steel X18H10T, the structural components of road-building, agricultural machines and ships. Accord- 
ing to the scientific research results in 1985, the monograph “Equipment for pulsed arc consumable electrode welding” 
was published in Energoatomizdat. 

The developments are protected by 17 copyright certificates and patents for inventions. Based on the research 
results, 31 papers were published in journals included in the list of the State Commission for Academic Degrees and 


Titles of the Russian Federation. 
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Brake rigging dynamic simulation under braking on a track section with ry 


irregularities (the case of a passenger car) 
I. A. Yaitskov, V. V. Kosarevskii 


Rostov State Transport University (Rostov-on-Don, Russian Federation) 






Introduction. The paper considers simulation of the dynamic processes of the brake rigging of a passenger car under 
braking on a track section with irregularities. The work objectives include the development of a “rigging brakeblock - 
wheel working surface” contact module in a full-scale computer model of a passenger car in the “Universal Mecha- 
nism” software package; and a computer simulation of the braking operating mode from 50 to 32 km/h considering ver- 
tical and horizontal track rail irregularities for determining the mechanism of variation of the longitudinal acceleration 
of the brakeblock and its angular acceleration. The subject of the study is the force interaction of the elements and dy- 
namic processes in the brake system of passenger cars. 

Materials and Methods. A new “rigging brakeblock - wheel working surface” contact module, which provides the de- 
termination of the longitudinal and angular accelerations of the brake rigging of a passenger car, is proposed to the 
“Universal Mechanism” software package,. The simulated modeling of the brake rigging system of a passenger car with 
KVZ-TsNII type II trolleys equipped with shoe brakes is carried out. 

Results. A full-scale computer model of a passenger car, which includes the designed contact module “linkage brake 
pad - wheel working surface”, has been developed in the “Universal Mechanism” software package. The car is present- 
ed as a system of solids connected by elastic and dissipative elements. Using computer simulation, the operating mode 
of braking was reproduced under reducing the speed of a passenger car from 50 to 32 km/h considering vertical and 
horizontal irregularities of a railway track. The simulation result was the laws of change in the longitudinal acceleration 
of the brakeblock and its angular acceleration under braking in the above speed range. Their spectra of longitudinal an- 
gular acceleration of the brakeblock were constructed. It was determined that the presence of track irregularities affects 
the spectral composition of the accelerations. In addition, under the superposition of the bogie-frame pitching and 
bouncing oscillations, when moving along an uneven track, the rigging block can move up and down the wheel-working 
surface within a range of up to 50 mm. The simulation functionality of the dynamic processes of the brake system of a 
passenger car was expanded in the “UM-Loko” software package. 

Discussion and Conclusions. The results obtained can be used in the design of new rigging brake systems of passenger 
cars and modernization of existing ones at the engineering enterprises and railway-car repair works. This, in turn, 


should ensure uniform distribution of efforts across all brake rigging brakeblocks of a passenger car. 
Keywords: rigging, brake system, wagon, block, dynamics, simulation, braking, track with irregularities. 
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Introduction. The problem of efficient and reliable operation of the brake systems of passenger cars becomes 
even more urgent when high-speed traffic is introduced on the Russian railways. The main part of modern brake sys- 
tems is the lever transmission system of a passenger car, whose efficiency and reliability depend directly on the quality 
of design, operation, maintenance and repair of the brake equipment. The mechanical part of the brake system combines 


a rigging brake transmission, an automatic rigging transmission regulator and friction brake components (brakeblocks 
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and pads). One of the key requirements for the brake rigging system is to provide uniform force distribution across all 
brakeblocks. However, under the operating conditions, there is some variation in the efforts of pressing the brakeblocks 
on the wheelsets both within the wagon and within each bogie. The nonuniformity of pressing brakeblocks may be one 
of the basic reasons for their uneven wear. The subject of the study is the force interaction of elements and dynamic 
processes in the brake rigging system of passenger cars. The level of theoretical studies on dynamic processes of struc- 
tural elements of the traction rolling stock is quite high [1-12]. 

Materials and Methods. To clarify and confirm the results of theoretical studies, it is required to carry out 
computer simulation of the dynamic processes of the brake rigging system of a passenger car, which occurs under brak- 
ing on a flat section of the track. Major causes of the dynamic processes in the contact “rigging brakeblock — wheel 
working surface” may be fluctuations of the car underframe when moving along a track with irregularities. This is be- 
cause the brakeblocks and the brake rigging system components are structurally connected with the bogie frame, which, 
due to deformation of the axlebox suspension, moves with respect to the wheel pairs rolling along the rails. It is re- 


quired to consider the simulation of the braking process in the presence of uneven tracks. 














Fig. 1. An example of the vertical track rail irregularities 


Files of vertical and horizontal irregularities of track rails, which are a generalization of the measurement re- 
sults, are available in the UM-Loko software package. Fig. 1 shows vertical irregularities on a track of 1000 m long as 
an example. 

Research Results. The simulation results of the dynamic processes of the lever transmission of the brake sys- 
tem on a flat track are presented in Fig. 2. The picture of the “wheel — rail” contact patches for all four wheelsets of the 
car, in the presence of track irregularities, is shown in Fig. 3. Comparison of simulation results with Fig. 2 provides the 
conclusion that the presence of irregularities causes the emergence of vertical oscillations of the car underframe and, as 


a result, a significant change in the forces within the “wheel — rail” contact patch. 
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Fig. 2. “Wheel — rail” contact patches and force distribution (on level) 


When moving along a track with irregularities, the oscillations of the truck-frame bouncing have an amplitude 
reaching 18 mm; frequencies of 0.8 are visible in the spectrum of vibrations; 0.95 and 1.07 Hz. The last frequency is 


close to the natural frequency of the car-body pitching. 
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Pitching oscillations of the bogie frame have amplitude up to 0.0025 rad. In the spectrum, there are frequencies 


of 0.25; 0.8; 0.95; 3.7 and 3.9 Hz. 
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Fig. 3. “Wheel — rail” contact patches and force distribution (uneven track) 
Fig. 4 and 5 show a graph of the longitudinal acceleration of the rigging brakeblock under braking in the pres- 


ence of track irregularities and its spectral composition. 
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Fig. 4. Longitudinal acceleration of the rigging brakeblock under 
braking (track irregularities) 


Fig. 6 and 7 show a graph of the angular acceleration of the rigging brakeblock under braking in the presence 


of track irregularities and its spectral composition. 
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Fig. 5. Spectral composition of longitudinal acceleration of rigging brakeblock (track irregularities) 
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Fig. 6. Angular acceleration of rigging brakeblock (track irregularities) 
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Fig. 7. Spectral composition of angular acceleration (track irregularities) 


The simulation results show that in the presence of track irregularities, the longitudinal and angular vibrations 
of the brakeblocks occur in the frequency range up to 25 Hz, and their spectral composition is somewhat different from 
the case when braking occurs on level. 
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Discussion and Conclusions. A full-size computer model of a passenger car is developed in the “Universal 
Mechanism” software package. A car is presented as a system of solids connected by elastic and dissipative elements. 
The structure of the model includes the developed “rigging brakeblock — wheel working surface” contact. Using the 
computer model, the operating mode of braking from 50 to 32 km/h is reproduced. The option when the track has verti- 
cal and horizontal irregularities of lengths of rails is considered. Because of the simulation, the variation pattern of the 
longitudinal acceleration of the brakeblock and its angular acceleration under braking from 50 to 32 km/h were ob- 
tained, and their spectra were constructed. The track irregularities affect the spectral composition of accelerations. In 
addition, under the superposition of bouncing and pitching oscillations of the bogie frame when moving along the une- 
ven track, the rigging brakeblock can move up and down along the wheel working surface with a swing up to 50 mm. 
The results obtained can be used under the design of new and modernization of existing brake rigging systems of pas- 
senger cars at machine-building enterprises and car repair enterprises to provide an even force distribution across all 
brakeblocks. 
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Introduction. Since welding is the only means to connect pipe lengths into a continuous line when constructing main 
pipelines, modern quality management systems for the welding industry products are based on minimizing the 
occurrence of specific defects. This is achieved through monitoring and documenting welding procedures. 

Materials and Methods. The analysis of monitoring systems customized for manual, mechanized and automatic orbital 
welding has shown that the industry urgently needs systems that not only control and document the welding process, but 
also predict the quality of weld joints. This actualizes the need to develop an intelligent module that could, basing on 
real-time monitoring results, predict the quality of welded joints on the fly. 

Results. Since the theoretical connection between the forecasting results and weld quality attributes is characterized by 
the interaction of a significant number of physical phenomena continuous in time, the results of welding can be 
described only by a sufficiently complete nonstationary physicomathematical model of the welding process. However, 
in order to be able to predict the results of welding directly during the monitoring of the process, a simplified 
forecasting model is proposed whose key feature is the ability to perform calculations synchronously with the real 
process, which is implemented in a real-time mode with a given interval. 

Discussion and Conclusions. The major obstacle to the successful functioning of the operational forecasting module, 
apart from the length of the numerical solution of equations, is an estimation error. To ensure the minimum error of 
virtual display during simplification, it is necessary to conduct comprehensive studies of the significance and influence 
of individual factors and phenomena on quality attributes. These observations determined the content and sequence of 
work on the creation and implementation of an intelligent module for the operational forecasting of welding quality. 
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also predict the quality of weld joints. This actualizes the need to develop an intelligent module that could, basing on 
real-time monitoring results, predict the quality of welded joints on the fly. 

Results. Since the theoretical connection between the forecasting results and weld quality attributes is characterized by 
the interaction of a significant number of physical phenomena continuous in time, the results of welding can be 
described only by a sufficiently complete nonstationary physicomathematical model of the welding process. However, 
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forecasting model is proposed whose key feature is the ability to perform calculations synchronously with the real 
process, which is implemented in a real-time mode with a given interval. 
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apart from the length of the numerical solution of equations, is an estimation error. To ensure the minimum error of 
virtual display during simplification, it is necessary to conduct comprehensive studies of the significance and influence 
of individual factors and phenomena on quality attributes. These observations determined the content and sequence of 
work on the creation and implementation of an intelligent module for the operational forecasting of welding quality. 
Undoubtedly, the information on the forecasting of the weld joint quality should enter a higher-level pipeline quality 
management system, as well as be analyzed by construction organizations in order to develop preventive measures to 


improve the organization and performance of welding work. 
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Introduction. The consistent development of the oil and gas industry infrastructure requires continuous im- 
provement of pipeline systems since pipeline transport is currently most preferred for transportation of hydrocarbon raw 
materials and products of its processing. It is known that at the present stage of technological development of the con- 
struction of main pipeline systems, welding is the only way to connect individual pipes into a continuous line on site. 

Given the consequences of accidents and failures in operation, trunk pipelines are classified as hazardous pro- 
duction facilities. Since, according to Rostekhnodzor, over 85% of accidents and catastrophes at facilities occur due to 
depressurization or failure of welded joints, the quality maintenance of welded joints is quite a chellenge. Therefore, the 
quality of welding is the basis for the safe operation of any pipeline transport system [1]. 

The quality of the weld is assessed by the dimensions of its cross section, the mechanical properties of the weld 
metal and the heat-affected zone (HAZ), the distribution of stresses and residual deformations, probability of cold and 
hot cracks, the presence of pores, non-fusion and other defects. Modern quality management systems for welding prod- 
ucts are based on minimizing the probability of specific defects. For this, accident-preventive measures are taken to 
prepare and implement welding processes. This approach contributes not only to improving the quality of welded joints, 
but also to improving welding technologies, rational selecting welding materials, and developing methods for monitor- 
ing welding processes. In this regard, we consider in more detail the basic welding technologies used in the construction 
of pipelines. 

Materials and Methods. Currently, the most commonly used method of joining pipes of main pipelines into a 
string is fusion arc welding. In the construction of pipelines, manual, mechanized and automatic welding is used. Manu- 
al arc welding is characterized by simplicity of implementation, equipment mobility, but it is quite laborious and re- 
quires a large number of qualified personnel. In addition, under manual arc welding, a significant number of defects can 
occur. Mechanized (semi-automatic) pipe welding, in comparison to manual arc welding, is more productive. However, 
it is not free from shortcomings, the major of which are increased spatter of electrode metal, problems of gas protection, 
aerosol emission, especially under FCAW wire welding. It should be noted that the methods of manual and mechanized 
welding are characterized by a high degree of subjective influence of performers. The “human factor” is less significant 
for automatic welding methods. However, orbital welding, in comparison to manual and mechanized welding, is less 
mobile. Currently, among installations for orbital welding, it is required to single out equipment for consumable elec- 
trode welding with controlled droplet transfer of electrode metal [2]. To reduce the impact of characteristic disturb- 
ances, it is very promising to use adaptive technologies that promptly correct the welding process [3]. 

Among consumables widely used for automatic and mechanized welding, gas-proof flux-cored wires can be 
distinguished. FCAW wires, despite a number of advantages, are used in much smaller volumes. 

In recent years, attempts have been intensified to introduce more efficient automatic welding technologies un- 
der construction, for example, plasma, laser, butt-resistance welding, and also welding with a combination of various 
heat sources [4]. However, all these welding methods are not yet industrially usable or are used on a very limited basis. 
The introduction of automatic welding processes reduces the human impact on the quality of weld joints. However, 
close attention is paid to issues related to the fulfillment of the prescribed requirements and instructions on the work 
execution by a welder or operator of welding equipment. 

Currently, welding monitoring procedures have become an integral part of scientific and technical support for 
the construction of trunk pipelines. Monitoring procedures provide observation of energy parameters of welding pro- 
cesses, prevention of causes for non-compliance of welded joints with the requirements of standards and technical doc- 
umentation (STD), and control of their elimination. At the same time, monitoring procedures provide the implementa- 
tion of a number of requirements for documenting the work performed. This is achieved by connecting special recorders 


with appropriate software to the modern welding equipment [5]. 
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A typical operation scheme of such systems is shown in Fig. 1. 
A number of domestic high-tech science-based enterprises are involved in the development of such systems of 
monitoring, registration and documentation of the parameters of the welding process. Systems can be integrated into 


modern digital sources of welding current or manufactured as separate units. 
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Fig. 1. Information flows of welding process monitoring 


When monitoring the process of arc welding for recording parameters, the authors of [6] developed a polyweld 
(up to 64 welding points) recording system for measuring, displaying and saving data on the welding current, arc volt- 
age and temperature of the weldable workpieces being welded. The registers of such a system can be located at a dis- 
tance of up to 300m from the installation site of a personal computer with dedicated software. 

Similar solutions are used by other developers of welding process documentation systems, for example, 
“Storm” enterprises, NPF ITS, “Telma’”, “Alloy”, and others. The “Weld Telecom” monitoring and control system de- 
veloped by “Alloy” provides data both from a single point and from a whole stock of welding equipment through a re- 
mote server using a wired or wireless Wi-Fi network. Such capabilities of the Weld Telecom system provide monitoring 
of the operations performed by the welders, and monitoring of the technical condition of the welding equipment. More- 
over, the system allows both monitoring the process and transmitting commands for adjusting welding modes, which 
should be considered a significant step in implementing the Industry 4.0 concept using sensors and networks. In addi- 
tion, the system provides documentation of welding processes with the automation of the collection, systematization 
and data storage for the formation of the weldable product certificates and their application in managing product quali- 
ty. The Weld Telecom system capabilities for processing and visualizing the welding mode parameters are shown in 
Fig. 2. 





Fig. 2. Weld Telecom system visualization options 


Foreign developers and manufacturers of welding equipment are actively involved in the development of sys- 
tems for registering and documenting welding processes. Fronius (Austria), Kemppi (Finland), Miller Electric (USA) 
and a number of other foreign companies demonstrate solutions similar to Russian companies for registering welding 


process parameters. 
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It should be noted that Production Monitoring systems from Lincoln Electric (USA) and Merkle Quality Con- 
trol from Merkle Group Inc (USA) have somewhat greater capabilities. These systems not only monitor and record the 
welding process parameters, but also identify and record data on their deviations [7]. The on-line Production Monitor- 
ing system compares the voltage, current, electrode wire feed speed, time interval of the welding process predefined by 
the user to the actual values. Merkle Quality Control system provides documentation and control of up to 8 parameters 
of the mechanized welding process with the possibility of archiving. Such capabilities of the deviation fixation system 
are provided by independent sensors of current, voltage, wire and gas. However, such systems work stably only with the 
proprietary equipment. 

Therefore, it can be stated that the most popular systems for monitoring welding processes are used to docu- 
ment their parameters, as well as to analyze the consumption of materials, use uptime, register the out-of-tolerance con- 
dition, and discipline welders. However, none of the systems answers the key question — what consequences will the 
identified deviations lead to. Given the entry, we need systems that not only control and document the welding process, 
but also predict its results. Based on the foregoing, we need systems that could not only control and document the weld- 
ing process, but also predict its results. Despite the attempts to predict the weld joints operational quality according to 
the monitoring results of the welding process parameters undertaken both earlier this century [8] and later [9], no fore- 
casting systems on site exist to this date. All this actualizes the need to develop an intelligent module that can quickly 
predict the quality of weld joints. However, such a module can only be created through linking the quality indicators of 
welded joints with the actual parameters of the welding mode recorded in the on-line mode. 

Discussion and Conclusions. Unfortunately, the theoretical relationship between the forecasting results and 
weld quality indicators is characterized by the interaction of a significant number of physical phenomena that are con- 
tinuous in time since they determine the heat supply, the conditions for the formation and crystallization of the weld 
pool, the dimensions of the weld and the HAZ microstructure [10]. However, the representation of welding processes in 
the form of complex, multifactor systems provides using mathematical modeling under their study. It is easier to collect 
the data required for forecasting under stationary shop-floor conditions [11]. However, under the conditions of the 
route, it is difficult to measure a number of the manual and mechanized welding parameters. For example, movements 
of a welder's hand and parameters of the electrode oscillations in cutting under manual or mechanized welding are un- 
controllable values. They can be judged only by indirect signs. 

The impact of the joint assembly on the quality of the weld formation should be considered, since even the as- 
sembled joints accepted by the Quality Control Dept (QCD) will have deviations within tolerance. If under orbital weld- 
ing, scanning laser-television systems can be used to determine the real profile of welded edges [12], then it is rather 
difficult to use them for manual and mechanized welding. Therefore, predicting the quality of manual and mechanized 
welding only by the results of processing the actual values of the process energy parameters can only be approximate, 
just a matter of judgment. To increase the reliability of such a forecast, real data on the distribution of grooving sizes 
along the pipe joint are required. In cases where it is impossible to use laser-television scanning, other approaches 
should be used. Thus, if we neglect the change in the gap under welding due to thermal expansion and deformation of 
the metal, then monitoring the joint assembly under welding can be replaced by measuring its dimensions before weld- 
ing starts. Since the gap size varies relatively slowly along the joint, it is sufficient to measure at several points, and 
determine the remaining dimensions of the joint through interpolating the available results. There is another possibility 
of increasing the reliability of the forecast when it is difficult to consider welding parameters. Their possible spread (for 
example, wire diameter tolerance, possible variations in welding and wire feed speeds, joint clearance, etc.) can be tak- 
en into account. However, in this case, it is rather difficult to impose deviations from the instability of those parameters 
that are not measured on the simulation result. 

It should be noted that due to fluctuations in the welding process parameters, the quality indicators of welds are 
unevenly distributed both along the length and thickness of the welded joint. Therefore, the final task of the module for 
predicting the quality of joints is to identify precisely those sections of the weld where the deterioration probability of 
quality indicators is critically high. Obviously, such a complex relationship of quality indicators and the results obtained 


requires using modern forecasting tools. 
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With this in mind, the predicted results can only be described by a dynamic (non-stationary) physicomathemat- 
ical model. Considering the problems of fixing an extended set of parameters, surrogate optimization techniques should 
be used in the model [13]. Given these considerations, Fig. 3 shows the structure of the dynamic physical and mathe- 
matical model used as part of the operational forecasting module for the quality of weld joints by the minimum number 


of analyzed parameters. 
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Fig. 3. Structure of physicomathematical model for predicting welding process results as part of intelligent module 


Before welding, information on the process is input into the program: the type of connection and the form of 
cutting, the thickness of the weldable parts, the grade of steel, the grade and diameter of the electrode wire, shielding 
gas, and the recommended welding modes. Then, this information is given in the certificates for welded joints, which 
also contains data on the joint number, the serial number of the welding equipment, the name of the welder, the start 
and end times of the joint welding process. 

A unique feature of the model for performing a quality forecast under welding is the need for synchronization 
with the real process of calculations, which is implemented in a real-time cycle with a given step determined by the 
weld pool response time. Therefore, it should be less than the length of the change in the depth of the crater on the sur- 
face of the bath when the arc current changes. The modeling of the thermodynamic state of the joint should be carried 
out at each step of the real time cycle since a major feature of arc welding, in addition to the possibility of the arc pene- 
tration into the crater of the weld pool, is its volume variation. With this in mind, the energy and mass balance should be 
achieved in a time not exceeding the selected step. The thermodynamic state of the metal makes it easy to determine the 
size of the weld pool, the thermal welding cycle, and the chemical composition of the weld, which enables to calculate 
the number of structural components and evaluate the mechanical properties of the weld metal and HAZ. In addition, 
using the known methods [14], stresses and strains in the vicinity of the weld pool can be calculated. 

Such solutions provide to evaluate not only the possibility of hot cracking, but also the tendency of the weld 


metal and HAZ to cold cracking according to the results of calculating the amount of martensite. 
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Embedding a physicomathematical model in the computer program of the operational forecasting module al- 
lows you to evaluate the probability of defects directly under welding (Fig. 4). 

For reliability of the results, it is required for the speed of virtual reproduction of the process to be greater than 
its real course. Therefore, the major obstacle to the successful functioning of the operational forecasting module is the 
duration of the numerical solution to the equations of the physicomathematical model. This circumstance forces us to 
simplify both the model itself and its numerical implementation. The model simplification can be achieved through lim- 
iting the scope of its application (specialization), as well as through reducing the number of measured parameters and 
determined quality indicators. However, to reduce the error of virtual reproduction under such simplification, it is nec- 
essary to conduct comprehensive studies to assess impact of individual factors and phenomena on the quality indicators 


of welded joints. 
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Fig. 4. Structure of system for predicting welding results: 

X are set parameters of welding process; S are controlled parameters of welding process; s are measurement results of 
welding parameters; M is data on joint geometry, cutting quality and spatial position; m are measurement results; 
Mp is data on physical properties of welded metal; Wq is regulatory data on quality requirements for weld; 

Sp is comparison base of mode parameters; K is information flow; 

R is signal of results of comparison of weld quality to regulatory requirements 


It stands to reason that the structure of the intelligent module for predicting the welding process results, its 
software will be specified and revised according to the results of pilot industrial use under the real-time forecasting of 
the quality of welded joints directly during welding work. 

The desired sequence of work on the creation and implementation of an intelligent module for real-time fore- 


casting of welding quality and the sequence of their implementation are presented in Fig. 5. 
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1. To develop the concept, structure and desired software of the 
intelligent module for real-time forecasting of weld joint quality 


2. To create a sufficiently complete dynamic physicomathematical 
model of the formation of welds under pipe welding 


3. To check adequacy of the physico-mathematical model according 
to the test samples welding results and make the desired updating 





4. To study impact of the welding mode parameters on the formation 
of defects and summarize the results in the form of simple 
mathematical relationships and dependencies 






5. To simplify the computer model to commensuration of the speed of 
virtual reproduction of the welding process with the real process 
reaction rate 








6. To validate and establish rules (algorithms) for decision making 
based on the results of real-time forecasting of weld joint quality 





Fig. 5. Content and operating sequence on creation and implementation of an intelligent module 
for real-time forecasting of welding quality 


To diagnose the technical condition of the main pipelines, it is required to collect all available information 
about the facility. To do this, information on predicting the quality of welded joints should go to a higher-level pipeline 
quality management system, for example, to an integrated diagnostic monitoring system for the linear part of main gas 
pipelines [15], or to a system for information analysis on the state of oil pipeline components [16, 17]. Construction 
organizations need information not only about defects in welded joints, but also about the causes of their occurrence. 
Therefore, data on the real-time forecasting of the quality of welded joints should be analyzed by construction organiza- 
tions to develop preventive measures to improve the organization and performance of welding operations. It is advisable 
to carry out such an analysis under stationary conditions, with the inclusion of a more complete computer model of the 
welding process and an information storage device with appropriate filters in the structure of the forecasting module. 
However, this approach also requires additional study on the significance of individual factors and phenomena by quali- 
ty indicators. Only then will the information on real-time forecasting the quality of welded joints become a truly effi- 


cient tool for taking preventive measures and removing the causes of defects. 


Conclusions 

1. A method is proposed for processing data for monitoring the welding process using a deterministic physi- 
comathematical model that provides a sufficiently accurate prediction of the welded joints quality directly during weld- 
ing, based on the relationships between the mode parameters and specified quality indicators of the joints. 

2. Real-time forecasting is provided by a special intelligent module whose software includes a computer pro- 
gram for the implementation of a physicomathematical model of online forecasting of the quality of welded joints. 

3. Since the major obstacle to the application of the proposed method of online data processing is high re- 
quirements for the speed of solving the model equations, then, to provide a high speed of virtual reproduction of the 
welding process, it is proposed to conduct comprehensive studies of the significance of individual factors and phenom- 
ena on quality indicators. These considerations determined the content and sequence of work on the creation and im- 


plementation of an intelligent module for the online forecasting of welding quality. 
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Introduction. The study of the austenite magnetic state in steels has provided the mechanism of the external magnetic 
field impact on steel under the hardening process. Previous studies have established a positive practical effect of heat 
treatment in a magnetic field. The work objectives were to create a computer model of the magnetic state of carbon steel 
austenite; to conduct computational experiments with a system of spins at various values of temperature and external 
magnetic field. 

Materials and Methods. The positions of the Ising model were used. The canonical ensemble of spins was modeled by 
the Monte Carlo method using the Metropolis algorithm. 

Results. The algorithm was implemented with the initial parameters selected through experimental data on the magnetic 
state of austenite. The inhomogeneity of this state without exposure to a magnetic field was studied. Data on the sizes of 
ferromagnetic clusters in austenite at various temperatures were obtained. It has been noted that the presence of an ex- 
ternal magnetic field counteracts the temperature disordering of spins. Data on an increase in the size of ferromagnetic 
clusters under growing magnetic field strength were obtained. 

Discussion and Conclusions. A two-dimensional computer model of the spin state of austenite of carbon steel has been devel- 
oped. The computational experiments with various parameters of the model have shown that there is a short-range order in the 
arrangement of spins above the Curie temperature. With a rise of the temperature of the system, the sizes of ordered regions 


decrease; but when an external magnetic field is applied, they increase. 
Keywords: quenching, steel, magnetic field, martensite, martensitic transformation, austenite. 
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Introduction. The magnetic state of austenite in steels and the effect that an external constant magnetic field 
can exert on it is considered. Previous studies have shown that heat treatment in a magnetic field has positive practical 
value [1-4]; therefore, it appears important to study the mechanisms of such an effect. C. Zener [5] pointed out the need 
to take into account the role of magnetic transformations when considering martensitic transitions in iron-based alloys. 
In calculating the thermodynamic characteristics of the martensitic transformation, he used two parameters — the mag- 
netic and non-magnetic components of the change in the free energy of pure iron. 

A number of researchers [6—8] found a deviation of the temperature dependence of the inverse magnetic sus- 
ceptibility from the linear Curie — Weiss law in the paramagnetic region of some steels. In this case, the curve was satis- 
factorily described by the Langevin function for superparamagnetics. This change in the magnetic properties of austen- 
ite in steels is associated with spontaneous fluctuations of the long-range ferromagnetic order. 

The experiments [9, 10] proved the existence of short-range order of spins (the so-called “swarms”’) above the 


Curie temperature using magnetic neutron diffraction. 
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A constant magnetic field during quenching cooling affects regions with an ordered arrangement of spins in 
austenite and their interaction. As a result, the nucleation of martensitic crystals on ferromagnetic clusters is initiated. 

Materials and Methods. At the moment, there is only one accurate and detailed method for experimental studies of 
“swarms”? of spins — magnetic neutron diffraction. In this paper, the magnetic state of austenite of carbon steel is studied through a 
computational experiment. The positions of the Ising model [11—13] are applied. Using the Monte Carlo method [14-16], a canoni- 
cal ensemble’ of spins was modeled, the use of which provides simulating the behavior of the spin system at a constant temperature. 
To obtain an arbitrary non-uniform probability distribution, the Metropolis algorithm is used [17] — a special case of the sampling 
procedure by significance when some possible samples are dropped. 

The Ising model is used to simulate phase transitions in magnetic substances or binary alloys. This lattice mod- 
el takes into account interactions between the nearest nodes. Spins are represented by the magnetic moments of atoms in 
the lattice sites, which interact with each other and with an external magnetic field (if any). 

The Ising model is based on the following simplifications: 

- kinetic energy of lattice sites is not considered; 

- when calculating the energy of interaction of spins, only the nearest neighbors are taken into account; 

- only two possible spin states are provided (positive f or negative | direction along z axis). 

As shown in [18], the study of the classical two-dimensional Ising model reveals common behavior patterns of 
magnetic systems near the phase transition temperature, even despite the simplifications. 

Using the Ising model, such macroscopic system characteristics as the average energy(E), average 
magnetization(M), specific heat C and magnetic susceptibility y were studied. When calculating the average values for 
all system configurations, the period until the system reached a relaxation state, which was excluded from the 
calculations, was taken into account’. 

Research Results. Let us describe the total energy of the system of two-dimensional spins in the Ising model 
taking into account the possible magnetic field h: 

E=—-JO jy Sisy — ADL Si. 
Here, J is the exchange coupling constant, which characterizes the strength of interaction of neighboring spins; N is the 
number of all spins; (i, 7) means summation over all the nearest pairs of spins. The number s; is associated with each i- 
th lattice site. It characterizes the direction of the magnetic moment and can take the values either +1 (if the spin is ori- 
ented in the positive direction of the z axis), or —1 (if the spin is oriented in the negative direction of the z axis). 

If the value of the exchange interaction constant is greater than zero, then the unidirectional state of two neigh- 
boring spins is energetically more favorable, i.e., the state with the lowest total energy is ferromagnetic. Otherwise, it 
will be more preferable that the neighboring spins are antiparallel with each other (antiferromagnetic state). 

The application of an external magnetic field along the z axis adds or subtracts additional internal energy of the 
spins according to their direction along the z axis. 

In further discussions on the thermodynamic characteristics of the system, the energies J and h will be meas- 
ured in temperature units. This is convenient when considering the interaction between spins, because when heated, the 
communication system between them weakens. 


We find the relation between the specific heat C and the system energy fluctuation in the canonical ensemble: 


C =5((e)? - (E?)). 


' The canonical ensemble is a statistical ensemble corresponding to a special physical system. It exchanges energy with the environment (thermostat) 
being in thermal balance with it, but does not exchange matter since it is separated from the thermostat by a membrane that is impermeable to parti- 
cles. For a brief description of such a system, two parameters are used: number of particles N and average energy E [13]. 

? Further in the text, when it comes to system configurations, the configurations after the system reaches relaxation state are implied. 
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The magnetization M of the system was calculated by adding all the values of s,; in this configuration. Then, 
the average magnetization (M) over all system configurations was calculated 

Magnetic susceptibility y at a given temperature: 

Xx = 5 ((M?) — (M)*). 

Possible configurations of the spin system will be determined by the values of 2” spin numbers s. Using the 
Metropolis algorithm, spin configurations S with probability w(S) can be generated, and then the sought quantities over 
all configurations can be averaged. The configurations that differ among themselves by one spin flip were considered. 
The decision on the flip of one or another spin (1.e., on the adoption of a trial configuration Sp) depended on the ratio of 
weight functions: 

w(S) : 
where E, and Ey, are energies of the systems with spin configurations S and S,, respectively. 

Spin sj flipped over and a new configuration was adopted if r> 1 or r< 1, but greater than a random number 
generated due to the uniform distribution over the interval from 0 to 1. Otherwise, the spin remained unchanged. At one 
Monte Carlo step, a number of flip attempts equal to the number of spins of the system N is made. 

For the two-dimensional case, toroidal boundary conditions were chosen: the lattice is represented by a ring in 
which spins located on the right boundary of the original square lattice interact with spins located on the left boundary. 
A similar ring interaction is provided for the upper and lower boundaries. This provides the same number of interactions 


for all spins. The interaction of any spin with neighbors can be considered as the interaction with one spin whose value 











is equal to the sum of the values of four neighboring spins (it can be 0, + 2, or + 4). In the two-dimensional case, the 





minimum possible value of the energy change when the considered (central) spin is overturned is 4/. 

Algorithm and initial parameters for the implementation of the model of the magnetic state of austenite. 
The algorithm structure and the basic functions developed under its implementation are performed in a mathematical 
package taking into account the recommendations [19]. 

The model was implemented through sequential performing the following operations. 

1. Setting the initial conditions: the number of spins of the system N, the exchange interaction constant J, the 
intensity of the external magnetic field h, temperature 7, the number of steps of the Monte Carlo method WN, (analogue of 
time). 

2. Performing the function that creates the initial configuration of the system considering the given N, J, h. 

3. Calculation of the instantaneous system configurations for each step N, considering N, J, h and the initial 
state of the system. 

4. Visualization of: 

- instantaneous system configurations at points of interest, 

- the dependence of the instantaneous energy of the system on time, 

- the dependence of the instantaneous magnetization of the system on time. 

Through observing spin configurations, the current magnetic state of the system can be evaluated. A study of 
these time dependences provides estimating the required relaxation time of the system t for considering it in subsequent 
calculations. 

5. Calculation of the average energy of the system(E), he specific heatC, the magnetic susceptibility y and the 
average magnetization(M). 

6. Calculation of parameters characterizing the sizes of ferromagnetically ordered regions: average size, maxi- 
mum size, minimum size, root-mean-square size deviation. (The calculation is based on a set of spin configurations.) 

The algorithm was implemented at the initial parameters selected due to experimental data on the magnetic 


state of austenite [8, 20]. The temperature of the Curie point of iron-carbon austenite of U8 steel is ~ 180 K, 1.e., if it 
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were possible to keep the austenite lattice to this temperature, then below the Curie point, austenite would obtain ferro- 
magnetic properties. The exchange interaction constant of the two-dimensional model was selected so that the system 
experienced a magnetic phase transition at a given temperature. The value of the exchange interaction constant for the 


two-dimensional case does not coincide with the real value of this constant for the iron-carbon alloy. However, in the 








framework of the two-dimensional model, the following parameters were selected: N = 625, J=0.78, h =0, N, = 2500. 
Due to this, under modeling in the temperature range including the region of the Curie point of austenite, effects have 
been observed that are completely analogous to the effects of a real system experiencing a phase magnetic transition 
from the paramagnetic state to the ferromagnetic state, namely: 

— maxima of the specific heat and magnetic susceptibility at the Curie point, 

— sharp decrease in the average magnetization during the phase transition, 

— discrepancy between the magnetization curve and the specific heat and magnetic susceptibility curves up- 


on transition to the ferromagnetic state (Fig. 1). 
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Fig. 1. Behavior of steel U8 properties near the ferromagnetic Curie point of austenite according to the results 


of computational experiments 


Fig. 1 shows that the magnetization does not drop to zero at the Curie point. The remains of the “anomalies” 
above the Curie point were repeatedly observed experimentally [6—10, 19]. L. D. Landau [12] and other researchers [6, 
8, 9, 19] attributed this to fluctuations of order at T > ©. Fluctuations of spontaneous magnetization should be most 
pronounced near the Curie point, because they tend to infinity at the Curie point itself. 

Fig. 2 schematically shows the temperature dependence of the short-range order parameter o from [19]. 


0 0 Tt 
Fig. 2. Temperature dependence of the short-range magnetic order parameter [19] 


Near the Curie temperature, o decreases sharply, and at T > ©, it asymptotically tends to zero. From the point 
of view of short-range order, this “unsharpness” of the Curie temperature predetermines the following phenomenon: all 
anomalies of ferromagnetics have a steep maximum at the Curie point, however, at T > © they should not immediately 
go to zero, but only gradually decrease in accordance with the gradual destruction of order at close distances. 

As shown in [19], the exchange interaction energy often determines the forces that are noticeably acting pre- 
cisely at close distances. Therefore, the thermodynamically-equilibrium state of a ferromagnetic should be determined 
by the short-range order between spins. This is a special case of the general problem of the cooperative behavior of in- 


teracting atoms (for example, during a first-order phase transition). 
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Ferromagnetism is largely a quantum phenomenon [15]. Nevertheless, even within the framework of the clas- 
sical statistical approach, the theory considering order at close distances is certainly useful in illustrative and qualitative 
terms [15, 19]. 

L. S. Stilbans took into account the “magnetic order” at close distances in the Ising scheme [21] and theoreti- 
cally explained the “remnants” of the heat capacity anomaly at temperatures above the Curie point. The abnormal heat 
capacity is due to the necessity of destroying the order at short distances, which is partially conserved above the Curie 
point. 

This method was developed by S. V. Vonsovsky [22-24]. He theoretically showed the differences between the 
ferromagnetic and paramagnetic Curie points. The results of his work prove that the paramagnetic susceptibility near 
the Curie point is finite, and does not tend to infinity, as follows from the theory without considering short-range order. 

Such a great interest in the simplified Ising model is explained by the fact that it provides solving the funda- 
mental difficulties arising in the theoretical interpretation of second-order phase transitions. 

Analysis of regions with short-range ferromagnetic order in austenite using a computational experiment. 
In this section, we study the inhomogeneity of the austenite magnetic state without external exposure to a magnetic 
field. 

To visualize the system instantaneous spin configuration at selected moments, the spins are represented by 
white and black squares. White squares are spins oriented in the positive direction of the z axis (for example, under the 
impact of exchange forces and (or) an applied external magnetic field). Black squares are spins oriented in the negative 


direction of the z axis (Fig. 3). 
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Fig. 3. System instantaneous configurations at selected moments (7 = 100 K, 4 =0) 


Fig. 3 shows the system visualization at a temperature of 100 K (i.e., below the Curie point) without the action 
of a magnetic field at different moments. The following is an interpretation of these configurations. 

- At time moment 0, the system is in the initial state, which is initiated using a random value generator (fully 
disordered state). 

- At time 0, the system is in the initial state, which is initiated using a random value generator (fully disordered 
state). 

- The transition to moment | means that each spin of the system (N = 625) was considered to make a decision 
on its revolution according to the conditions described above. 

- By the third step, each spin was examined three times for a flip. 

- Subsequent moments show that the system tends to a state with a ferromagnetic order. 

- At the 30th step, a completely ordered state is recorded (with the exception of one spin). A similar picture 
(with small fluctuations of 1-3 spins) is observed at all subsequent moments. The total number of steps is N, = 2500. 
Consequently, the system has come to a stable state (that is, the relaxation stage has been completed before the 30th 


step). 
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It should be noted that the described picture is characteristic of the current generation of the model with given 
conditions. When generating the model, a random number generator is used, therefore, the results of the recounting will 
differ in details, but qualitatively, the model always behaves in a similar way under the same initial conditions. 

Consider the graphs of the instantaneous values of energy (Fig. 4 a) and magnetic moment (Fig. 4 5) of the 
system versus time. 


500 750 tag? 1.256105 1.50% 1.7510" 210% 2.25-10% 2.510% N, 








0 250 500 750 140? 1.2510 1.510% L750 240" 2.25108 2.56105 
b) 
Fig. 4. Dependence of the instantaneous values of the system total energy (a) 
and magnetic moment (4) on time 

Fig. 4 shows that the system quickly tends to an equilibrium state (relaxation process). After the transition of 
the system to the equilibrium state, only small fluctuations are noted. Upon reaching equilibrium for a given tempera- 
ture, macro parameters (average energy, magnetization, magnetic susceptibility, specific heat) were statistically calcu- 
lated from the system configurations. For this, the ratio of the relaxation time of the static system to the equilibrium 
state was estimated from the corresponding graphs. The initial system configurations were cut off for 400 steps (16% of 
the total calculation time), that, on the safe side, blocked the relaxation stage of the system, which did not exceed 150 
steps (see Fig. 4). 

When hardening U8 steel, the martensitic transformation starts at a temperature of ~ 500 K. It is interesting to 
study the model behavior near this temperature, because under quenching in a magnetic field, the field action initiating 
the phase transition is manifested here. Naturally, at this temperature, the paramagnetic state of austenite will be stable, 
as the model showed that. However, in the spin state of austenite, nanoscale regions with short-range order in the ar- 


rangement of spins are observed [9]; they play a significant role under the impact of an external magnetic field [8] (the 
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effect of the field is manifested in the influence on the size and stability of such regions, but they exist without an exter- 
nal field). 

To obtain the most reliable statistical data on ferromagnetically ordered regions in austenite, a system was sim- 
ulated with the maximum possible (according to the technical limitations of the RAM and processor speed) number 
(N = 10000) of spins at temperatures from 100 to 500 K. Other initial conditions remained the same. 

Fig. 5 presents the results of measuring ordered regions in austenite. The size was determined by the number of 
atoms per horizontal and vertical sections of the ordered region: L is average size, Max is maximum, Min is minimum, o 


is root-mean-square deviation. 
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Fig. 5. Temperature dependence of area sizes with ferromagnetic order in austenite 

According to the data obtained, fluctuations of the short-range magnetic order in U8 steel are preserved up to 
temperatures of about 500 K, when martensitic transformation starts under the steel quenching. In this case, the maxi- 
mum sizes of such regions in austenite at a temperature of 500 K can reach 20 atoms in the cross section (~ 7.2 nm). 
The maximum size at temperatures in the region of the Curie point and below is limited by the model conditions and 
does not exceed 100 atoms (i.e., the maximum transverse size of the modeled region). The observed oscillations on the 
curve of the maximum cluster size are associated with the use of a random number generator in the model. 

The Weiss theory considers only the order at long distances [19], which is determined by the difference in the 
number of spins oppositely directed throughout the domain, regardless of their detailed relative position. It turns out that 
the short-range order is simply equated to the long-range order, i.e., the numbers of the nearest neighbors of different 
orientations for any spin are assumed to be proportional to the total numbers of multidirectional spins in the domain. 
Thus, the local magnetization of the nearest neighbors at each spin is equal to the magnetization calculated over the 
entire crystal. A different picture is observed in a real crystal: parallel spins under the action of exchange forces tend to 
unite into separate “swarms” [9, 19] like gas atoms during its condensation or the formation of fluctuations in a non- 
ideal gas. As can be seen from Fig. 5, at low temperatures, these fluctuations are very large. At high temperatures, a 
more reasonable approximation does not consider the short range order. However, in this case, there are still fluctua- 
tions of parallel spins in small volumes caused by the exchange forces and leading to nonzero local magnetization in the 
absence of the resulting magnetic moment in the entire volume (i.e., in the absence of long-range order). 

Thus, the model predicts the existence of regions with short-range magnetic order in the austenite of U8 steel 
before the martensitic transformation. In these areas, with a high degree of probability, martensite nuclei arise. Clusters 
with ferromagnetic ordering can overlap defects in the crystal lattice and instability regions of the crystal lattice [4] that 
occur before the phase transition. In this case, the probability of martensite nucleation is especially high. 


Adding an external magnetic field to the austenite magnetic state model. The presence of an external mag- 
netic field counteracts the temperature disordering of spins. Fig. 6 shows the dependences of the cluster parameters in 
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austenite with short-range order (the parameters are the same as in Fig. 5) on the external magnetic field strength for U8 
steel at a temperature of 500 K. 
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Fig. 6. Dependence of sizes of regions with ferromagnetic order in austenite on strength of external magnetic 
field at a temperature of 500 K 

As can be seen from Fig. 6, with increasing magnetic field strength, cluster sizes grow up to the state of com- 
plete magnetic ordering at high magnetic field strength. Thus, under quenching in a magnetic field, regions with short- 
range order expand. Moreover, the probability of martensite nucleation is higher, and the number of nucleation sites is 
greater. 

Discussion and Conclusions. A two-dimensional computer model of the spin state of austenite of U8 carbon 
steel is developed. The computational experiments with various parameters of the model have shown that there is a 
short-range order in the arrangement of spins above the Curie temperature. With an increase in the temperature of the 
system, the sizes of ordered regions decrease; however, even at the temperature of the onset of martensitic transfor- 
mation (500 K), short-range fluctuations with a maximum size of ~ 7 nm are possible. When an external magnetic field 
is applied and its intensity is enhanced, the sizes of clusters with ferromagnetic ordering increase. Using the presented 
model in future research, it is supposed to study in more detail the impact of the magnetic field on the characteristics of 


the swarm of spins, their number, lifetime and interaction features. 
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Determination of linear characteristics of rotor mounting groups under load 
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Introduction. The paper considers analytical studies on the “rotor — gapped-type support” dynamic system under pro- 
cess loading. The research objective is to obtain expressions for determining the equivalent stiffness characteristics of 
the system. 

Materials and Methods. A rotor rotating in the elastic gapped-type supports is considered. A dynamic model that ena- 
bles to consider the problem of determining the linear equivalent stiffness characteristics of mounting groups is pro- 
posed. To solve the problem, a system of differential equations is compiled, and a detailed analysis is performed. 
Results. From the obtained dynamic equations of the system in question, we can calculate the static angular deviation of 
the rotor pins due to the action of the load. The proposed expressions for determining equivalent stiffness characteristics 
testify that it is possible to study the rotor dynamics as on the linear elastic supports with the above parameters. The 
obtained system of equations is analyzed, and all special cases of applying the first approximation formulas for equiva- 
lent stiffness of the rotor mounting groups are listed. 

Discussion and Conclusions. The results obtained make it possible to study many dynamic processes on the basis of 
linear differential equations considering the nonlinear properties of the system. For shavers used in the production of 
leather materials, the determination of rotor vibrations in the horizontal plane provides the quality and accuracy of oper- 


ations. 


Keywords: mechanical engineering, rotor, dynamic rotor equipment, rotor balance, rotor mounting groups, linear equiv- 
alent stiffness parameters of rotor units, rotor operating modes. 
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Introduction. Technological progress imposes new, higher requirements for the quality of engineering prod- 
ucts. Mainly, it is about high reliability, durability of machines, their productivity and safety. All these parameters 
should be taken into account and calculated at the design stage. 

From the point of view of production processes, the dynamic rotor equipment that provides the continuity of 
the process is of particular importance [1-3]. In [4], structural designs and the principle of operation of various rotary 
machines are presented; patent searches for each type are performed. In [5], the nature of the dependence of equivalent 
stiffness on the frequency of oscillations for a specific model is considered; the dependence of the critical frequencies 
on the equivalent stiffness of the supports is given. In [6], the authors consider in detail the issues of vibration control, 
vibro-adjustment work and prevention of increased vibration of rotor machines, and equilibrium of rotors. In addition, 
vibration sources are listed here and basic information from oscillation theory is provided. The problems of the dynam- 
ics of a rigid unbalanced rotor with four degrees of freedom are discussed in [7—9]. The study [10] is based on the as- 
sumption that bearing reactions are quasilinear forces with cubic nonlinearity. With this in mind, the effect of the radial 
clearance on the movement in space of a dynamically unbalanced rotor under the impact of internal friction forces is 


considered. In [11], the interrelation of transverse and torsional vibrations arising under the rotation of the centrifuge 
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rotor is shown. A linearized mathematical model of the rotor in elastic supports is developed due to the impact of trans- 
verse and torsional vibrations. 

The literature review suggests that the study on the oscillatory process and the corresponding parameters of 
linear mechanical systems is of interest. In this paper, the rotor, in particular, its gapped-type supports, is considered. 
The study objective is to determine the values under load of linear equivalent stiffness parameters of the named mount- 
ing groups. 

Materials and Methods. A dynamic model that provides for rotation in elastic gapped-type supports was used 


as a rotor (Fig. 1). 





Fig. 1. Dynamic model of a rotor that rotates in elastic gapped-type supports 


In Fig. 1, the following designations are adopted: 1 is rotor support; 2 is the contact boundary of the support 
with the rotor axle; 3 is the curve of relative motion of the rotor axle center; m is the rotor mass; B (a) is the radial clear- 
ance in the rotor bearings (represented by the angular displacement function of its journals from the vertical guides); C,, 
C, are the general stiffness of the casings of the mounting groups in the horizontal and vertical guides, respectively; yc: 
Zcr is the static shift of the center of rotor mass in the horizontal and vertical direction caused by deformations in the 
mounting groups; a is the dynamic tilt angle of the rotor journals; o,, is the angle of inclination from the equilibrium 
position caused by the process load on the rotor; y,, Z, are total dynamic displacements of the rotor mass center in the 
above directions; f= Fz / Fy is the relationship of the vertical and horizontal components of the process load. 

To determine the potential (/7) and kinetic (T ) energies of the system under study, we use the following 


equalities: 


= et Ly. + Vo, ~5(4+0,,, )Sin (a+ o,,, )+5(O,,, )SiN O.,, ie + 


+ : C, [z + Zon — 8(O + OL, )COS(O+ OL,,, ) + 5(OL,,, )COS OL, i 
—mg| z, + Zem +8( 4, )COSO,, |; 
1 ; . 1 
T=sm(>; +2))+5 40. 


Here, A is polar moment of inertia, @ is the angular velocity of a rotating rotor. 
It should be borne in mind that the acting forces are not potential. Moreover, the generalized forces to which 
the previously introduced coordinates y,, Z,, @ correspond, take the form: 


0, =F,, 0, =-F, 
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5 [d(a +O,, )COS(a+O,,,) +5’ (a+ o,,, )sin(a+o,,, )] + 


7 +f[8(a+a.,, )sin (+c, )—9'(a+c,,, ) Cos (c+ c,,, ) | 


Here, S'(a + on ) is a angle derivative of radial clearance; f =— is the relationship of vertical and horizontal com- 


K 

ponents of the process load. 

Under these conditions, equations expressing the dynamics of the system can be written as: 
my, +C,-y, -C,8(a+o,,, )sin(ato,,,)+C,5(a,,, )sina,,, = 0, 


mz,+C_-z,-C,8(a+0.,,, )cos(a+o,,,)+C_8(a,,, )cosa,,, = 0, 


cm 


mg [3(a +0,,,, )Sin(a+a,,,)—8'(a+O,,, )COS(a+ O,, )] - 


-C,[ y, -8(a+a,,, )sin(a+a,,, )+5(c,,, )sina,,, ] 


x[ 8'(a+a,,, )sin(a+a,,,)+5(a+c,,, )cos(a+,,, ) ]+ (1) 
+C, [ Z, —8(a+a,,, )cos(a+a,,, )+5(0.,,, )cosa,, |x 
x| S(a+a,,, )sin(a+a,,,)-8'(a+a,,, )cos(a+ 6, )] = 


[ (a +a, )cos(a+o.,,,)+8'(a+a,,, )sin(a+o,,,) [+ 
. +f [8(a +O,,,, )Sin(a+ o,,,)—8'(4+O,,, )COS(@+ O,,, )] 
Having examined in more detail the third equation of the system (1), we determine the angle of inclination a, 
considering a = 0. After transformations in the equations (1), we obtain a simpler system: 


my, +C? 


KB 


eV Coe tg Oy 
mzZ,+C,, . +C2 -y, =0. ®) 
An analysis of the equations (2) enables, in the first approximation of the values, to consider the elastic 
gapped-type support as a support with linear elastic characteristics (both in horizontal and vertical guides): 
Fy Cog Vn POs Eas Fa = ogg Fg Hoe Yi 


Here, 
C, ‘r+cos’a,, (0+ f-2z:2°-C, C, ‘r+c08 O,,, (0+ f-2-0°-C, 
Core = [me : ( = NN ee -[me ’ ( : I 
A A 
(38 ae Gos Me BOF 2) 
A 
r=8 (a,,,)+28'(a,,) —5(a,,)-5"(a 
0=8(a,,,)+8'(G,,,)18O,, 5 2=8(G,,, )18O,,, —5' (A, ) > 


A=mg-r+cos> G., (0+ f+z):(C, “Z +C, 0°). (3) 











5) 


cm J? 


Chas Coe 


oK6 ? Ke 


Research Results. In the system (3), the quantities C” 


OK6 ? 


express the equivalent stiffness characteris- 
tics of the gapped-type support. 

To obtain the dependences of the equivalent stiffness characteristics of a non-linear mechanical system, we use 
the expansion of trigonometric functions in the system (1) taking into account higher-order summands. To determine 
the required values, we indicate: 

y, =A, sin@t, z, =B,-cosa@,t, 
where B,, A, is the amplitude of the general oscillation of the rotor in the vertical and horizontal guides, respectively; 


@, is the frequency of oscillations. 


Thus, 
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C? va [ ®, (y,,z,)sino,tdt, 


KE 


ee (4) 
Zz , “ 
«<=—J]@ ,»Z, )COs @,tdt. 
KE B J Zz (y, ) 1 


Here, @, ( VyoZq) , @, ( Vos%q) are functions of y,, z, coordinates. These quantities can be determined through the 


system of differential equations (1). 


Solving the system of equations (4) by integration methods, we represent the quantities in the form: 


C, (mgr +C. COS G.,,, -b-z") C? cosa, -y-0-b°(C; -0°- AY +C) +z" - B?) 
Po =: 








Cc” = 
sal A 8A° , (5) 
a C. (mgr +C, cosa,,, -b-8°) C? COS Og, W+z-b° (C?-0°- A? +C? +z? - B?) 
= + : . 
oa A 8A° 


Here, 
b=cos’ o,,, (0+ f-z), 


be 5(a.,,,)+ 38' (02, )tgo.,,, —38" (a, )—8"(a.,, )tgou,,, > 
& = 8 (0,4, ALO yy — 35" (Og, )— 38" (Oey, EZ FE" (Oey) 


Analyzing (5), we can talk about the interdependence not only between the stiffness of the housings in the rotor 


bearings C, , C,, but also the amplitudes A, , B, of its general vibrations. 


It must be emphasized that stiffness C””, is expressed as the relationship of rotor movements in the horizontal 
Ip p p 


and vertical planes. In the case C?°, = 0 in the system (2), the equations will not be connected, which enables to consid- 


er all possible options. They are listed below. 
— In the absence of changes in the radial clearance over time and a,, =0, we have C;,, = C,. This option is 


typical for idle rotor operation. Moreover, the radial displacement of the rotor with respect to the support in the vertical 


guides is a very small amount, although of a higher order than the displacement in the horizontal guides. 


— If we accept z =0, then C;,, = C_. This situation is characteristic of the radial movement of the rotor rela- 


tive to the support surrounded by a point through which the horizontal tangent passes to the trajectory of the relative 
movement of the journal center. 


— When the value of the radial clearance in the bearings is zero, we have C>,, =C, andC>,, =C,. 


9KB 


= 90°, then C’ 


— If there are no changes in the radial clearance and «.,, Sxe = C, . This is possible if the rotor op- 


erates at under the process load with a selected gap, i.e., in horizontal rails, the radial displacement assumes small val- 


ues in comparison to the vertical direction. 


— When 0=0, C:,, =C,,. In this case, the radial displacements of the rotor take place in the vicinity of the 


point through which the vertical tangent passes to the line of relative motion of the journal center. 


— It is possible to accomplish the equality 0+ f-z=0. This is true in the case of a radial displacement of the 
rotor with respect to the support surrounded by such a point through which a normal is drawn which coincides with the 


line of action of the resulting forces F, and F,. In this case, C}., =C, andC>,, =C,. 


r 


CG 


Consider the characteristics of the rotor on elastic supports the rigidity of which is equal to C? * . Let us 


oK6 ? 


denote ¢ the eccentricity of the rotor, and @ its angular velocity. We represent the oscillation equations in the form of 


the system: 


KE 


(6) 


ee +C!-y, =m-é-o° sinot 


mi,+C.,+Z, =m-E-@ cosot 


KE 
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Solving the (6), we write the expressions for the oscillation amplitudes: 


2 2 
m:&-O B= atta (7) 


C2 -m-o 


KE 


A 


y 





== 7 
Co, Mm: 


Note that the summands C” 


OK6 ? 


Ce in the denominators depend on the amplitudes 4,, B_. We represent the 


graphically considered system (Fig. 2). 


@ (@2)°,(@,)” 








Fig. 2. Graphical representation of dynamic characteristics of absolute oscillations of the rotor in gapped-type supports 


The surfaces F,, P,,, A,, PB, (resonant) and C,, C, (skeleton) are constructed under the assumption that 
the amplitudes A, and B, vary independently based on (7) and the equalities C}., = m(@? y » C= m(o; y : 
The amplitudes of forced oscillations A,, B, correspond to the frequency of disturbing forces @,. Given this 


fact, to obtain expressions of the indicated amplitudes, we first construct the plane w” = @, . 


At the intersection lines of the constructed plane with the surfaces F,,, P,, and P 


lz? 


ae P,, there will be points 


with identical coordinates. These points are the desired amplitudes of the rotor oscillations. 


Planes C i , Cu" are the skeleton surfaces of a linear system. In this case, only the first summands are taken in 


the expressions for equivalent stiffness. They indicate that the frequencies of free vibrations of a linear system are inde- 
pendent of amplitudes. 


Resonant surfaces represent rotor oscillations in horizontal guides, and skeleton surfaces — in vertical guides. 
They are built in various coordinate octants. The skeleton surfaces are elliptical paraboloids in shape. w’, m; represent 


the frequencies of free oscillations of the system in Fig. 2. 

Discussion and Conclusions. Multiple points may be suitable for a single frequency. This means that several 
modes of oscillations including unstable ones are possible in the studied system. For the transition of a system from one 
stable mode of motion to another, some external impacts are necessary, which is characteristic of nonlinear systems. 

As an example of the application of the dependences obtained, we can cite the problem on forced oscillations 


of the rotor due to its static imbalance. 
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Dependences presented in the paper (4) indicate that equivalent angular stiffness properties are interdependent 


not only through the stiffness characteristics C,, C_ of the bodies, but also through the amplitudes 4,, B, of general 


rotor oscillations. 

As a result of the study on the linear characteristics of the mounting groups of the loaded rotor, the following 
results are obtained. 

1. A dynamic model of a rotor rotating in elastic gapped-type supports is proposed. This model provides inves- 
tigating the problem of determining the linear equivalent stiffness characteristics of mounting groups. 

2. The system of equations is analyzed, and the possibilities of applying the formulas are listed. In particular, 


they can be used to determine rotor oscillations in the horizontal plane on shavers used to produce leather blanks. 
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Introduction. The paper considers the mathematical analysis of the fractional composition of dust generated during the 
operation of the rail-cutting machine. It is established that the studied polydisperse material is well described by the 
one-parameter exponential distribution. At the same time, the lognormal particle size distribution, whose parameters are 
determined by mathematical programming methods, seems adequate for the purposes of calculating cyclones. The work 
objective was to develop mathematical methods for correct averaging of the size and mass parameters of dust under the 
solid metal machining. 
Materials and Methods. We studied the possibility of approximating experimental data by Rosin — Rammler 
distributions (classical, generalized three-parameter P(x, D, n, m), and simplified exponential P(x), in which n = 1). The 
corresponding results were compared to each other and to the data of approximation of the lognormal and double 
lognormal functions. These results indicate close approximation quality using the following model distributions P(x): 

e five-parameter double lognormal; 

e three-parameter type of Rosin-Rammler; 

e two-parameter classical Rosin - Rammler; 

e one-parameter exponential. 
Results. The primary physical analysis of cutting waste was carried out by the laboratory measuring complex Fritsch 
Analysette 22 Compact which uses the LALLS — low angle laser scattering method. The built-in software provides 
output of measurement results in primary graphic and digital forms. It was found that the simplest exponential 
distribution is best suited for a detailed analysis of the dust particle-size distribution based on the experimental data. 
This distribution enables reproduction of all the integral indicators provided by the instrumental measuring complex 
along with the graphical data. 
Discussion and Conclusions. The results obtained can be used to rationalize the local suction machine, and the 


mathematical models and algorithms can be used for the parametric analysis of any dust captured by cyclones. 
Keywords: metal cutting, dust, size distribution, statistics, mathematical programming 
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Introduction. Dust generated during machining (cutting, drilling, grinding, polishing, etc.) of hard metals and 
alloys can cause direct or indirect damage to human health and lead to environmental pollution. To eliminate these 
negative consequences, cyclones are widely used — devices of general and local air purification using aerodynamic 
capture of dust by inertia forces with its subsequent screening from the air stream into the accumulator [1]. The 
efficiency of cyclones is ensured by calculation, the purpose of which is to guarantee the capture of solid particles 
suspended in an air stream of a given category at minimal economic costs. Since particle capture is provided by 
competition between inertial and aerodynamic forces, the corresponding physical criteria form the basis for the 
calculation of cyclones [2]. The key components of these criteria are the dimensional and mass characteristics of dust 
particles. Due to the natural heterogeneity of industrial dust, these characteristics are of a statistical nature, which puts 
forward narrow requirements for the correctness of their averaging when calculating dust cleaning systems (cyclones). 
Thus, reliable scientific information on the dimensional and mass parameters of dust for various types of mechanical 
metalworking is a challenge since their practical application provides for a rational organization of dust collection 
through cyclones. 

This study objective is to develop mathematical methods for the correct averaging of size and mass parameters 
of dust under the machining of hard metals. The results obtained were used to rationalize the local exhaust machine; and 
the mathematical models and algorithms developed by the authors are applicable without significant restrictions for the 
parametric analysis of any dust to be captured by cyclones. 

The theoretical framework for the study. The dust particle is under the impact of gravity and inertial forces, 
which is proportional to the mass of the particle equal to p-x°, where p is the density of the forming material in [kg/m*]; 
x is the characteristic particle size in [um]. Therefore, knowledge of the dust inertial properties comes down to 
knowledge of its density and characteristic size. The magnitude of the aerodynamic force acting on such a particle from 
the flow side is proportional to the square of its characteristic size x’ and does not depend on density. The ratio of these 
forces appearing in the cyclone efficiency criterion is, respectively, proportional to the complex p-x. However, two 
important circumstances should be taken into account: firstly, dust particles differ significantly in size and, secondly, 
the shape of each dust particle is unique and far from the standards used (sphere, cube, etc.). These features put forward 
very narrow requirements for the procedure of double averaging of the characteristic size of a dust particle by its 
component size and shape. Obviously, the technique averaging the characteristic size of dust particles is targeted: in our 
case, it corresponds to the calculation of cyclone efficiency. 

For the first time, J. Sauter studied systematically the parametric averaging of polydisperse media [3, 4]. The 
key results of his work are as follows: 

1. For various purposes, specifically averaged medium particle sizes from an inhomogeneous population are 
important. Since the average size is a value expressed in fractions of a meter, the following method of dimensional 
averaging of an ensemble of particles with a distribution function F(x) and, accordingly, a probability density P(x) = 
dF(x)/dx is obvious: 

<x>y = [J PQ) x! dx/| P(x) x dx] . (1) 

2. The formula (1) implies that all dust particles are characterized by a single size, i.e., they have the shape of a 
ball. Therefore, the value <x>j is called the Sauter diameter and is designated as D, (most often, the Sauter diameter is 
understood as D3). 

3. If the shape of the particles is substantially irregular and is characterized by two or two parameters, the 
shape factor is also introduced into consideration. The physical meaning and practical application of various Sauter 


diameters are given in Table 1, and important information about the shape coefficient is in [5]. 
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Numerous studies on the dispersion of various media and materials allow us to state [6] that the solid particles 
obtained as a result of single crushing are distributed in size according to the two-parameter Rosin — Rammler law: 

F(x, D,n)=1-e@", P(x, D,n) = I/D - (x/DY"!- (2) 
where the value <x> = D - G(1 + 1/n) characterizes the average particle size, and n is the degree of dimensional 
heterogeneity of the ensemble (the smaller n, the higher the polydispersity of the powder). 

Under multiple grinding, the powders consist of particles whose sizes satisfy the two-parameter lognormal Gauss- 
Kolmogorov distribution [6]: 


aly lgx-lgD ; 
Ige 2 Igo 


—— eC 
V2n-Igo-x 


In the distribution (3), the parameter lg D corresponds to the conditional average particle size, and the 


P(x, D, 0) = (3) 


parameter lg o corresponds to the spread in real particle sizes around the conditional average. 

The principal advantage of the Gauss —- Kolmogorov model is the convenience of recalculating the values of D, 
from the linear Hatch — Choute relations [7], which links them to the quantities D and o. The form of these relations is 
such that at any two known quantities D,, and Dy, all other quantities can be calculated. 

It is important to note that the powder analyzed by the authors (rail cutting waste) is not necessarily described 
by the classical models given here. Firstly, the cutting process contains elements of both uniqueness (each contact of the 
abrasive wheel with the rail material is unique) and multiplicity (such interactions are extremely repetitive). Secondly, 
the shape of metallic filings is far from spherical. Finally, along with metal filings, cutting waste contains particles of a 
chipping abrasive. The content of the latter, due to the characteristics of the process and the requirements for the tool, 
can vary significantly. Thus, the study on the distribution of exhaust dust by particle size seems to be a practically 
important and scientifically significant task. 

Experimental data. The initial physical analysis of the cutting waste was carried out using the Fritsch 
Analysette 22 Compact laboratory measuring complex through the LALLS — low angle laser scattering method [8]. 
The embedded software implements measurement results in primary graphic and digital forms. The disadvantages of 
the software part of the complex are: lack of detailed information on the conversion algorithms of the measured 
quantities and on the output data nature, as well as the quantitative errors of the data displayed in the form of graphs. In 
particular, it may seem that the scale of the differential distribution function (probability density) is given on printouts 
with a several-fold error. A deeper examination allows us to conclude that in fact on this graph, the dependence of the 
quantity P(x,) dx, on x; is shown, and the partition of particles by size into groups of width dx, is not uniform. In 
addition, the integral characteristics of the particle Dj size distribution calculated by the program are not documented, 
which requires their verification for compliance with both the original graphic data and the classical models of Rosin — 
Rammler and Gauss — Kolmogorov. 

The noted circumstances, when evaluating the results of the analysis output by the built-in program, determine 
the implementation of additional measures: normalization of the initial curve of the differential distribution function, as 
well as checking the concordance of the integrated indices D, of both this distribution function and the basic model 
distributions. The implementation of these tasks requires high-quality digitization of the graphical results of fractional 
composition provided by the Fritsch Analysette 22 Compact, and this requires the specialized software [9, 10] and the 
development of appropriate data verification algorithms. 

Digitization technique and experimental data verification. To digitize the initial graphic data of the 
dispersion analysis obtained using the Fritsch Analysette 22 Compact instrument, the specialized Grafula program [11] 


was used. This information freeware provides for translating graphically represented dependences into a tabular form. 
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The digitization procedure is reduced to reading the graph, placing the axes of the Cartesian coordinates on it, setting 
the scale and drawing a sufficient number of markers on the line representing the dependence. 

The result of automatic digitization refers to the position of the points input by the user and is formed into an 
Excel-compatible table. The error in the resulting digital data has several components: 

* error in the schedule generation due to the features of the Fritsch Analysette 22 Compact complex; 

* defects in printing or displaying as an electronic photograph (screenshot); 

* inaccuracy in digitizing the graph through the Grafula program; 

* impossibility of accurately marking the curve on the graph due to limitations both in the resolution of the 
system and in the psychomotor capabilities of a person. 

The last of these errors, apparently, is the most significant. The noted circumstances require supplemental checks 
of the digitization results to eliminate critical errors and evaluate correctly the result error. 

To verify the digitization results, this study was carried out in several stages. At first, as applied to data in a 
normal particle size representation relative to those presented on a logarithmic scale, and then the initial analysis results 
were compared graphically. In this case, both digitized series were renormalized in order to eliminate the systematic 
error and ensure that the condition important for the distribution function F(x) is fulfilled: F(oo) = 1. Subsequently, on 
the basis of each of the results, the uniform model distributions were constructed whose empirical parameters were 
determined by mathematical programming methods [12] and compared. The last test stage was comparison of the 
averaged indicators of the dust dimensional composition calculated from the constructed model distributions with the 
integral characteristics of the distributions that are issued by the Fritsch Analysette 22 Compact instrument software. 

It turned out that the distortion of the distribution function F(x) as a result of the inability to extract it from the 
instrument directly in the digital representation is 15% when digitizing the graph in a linear (particle size) scale and 8% 
when digitizing the graph in a logarithmic scale. Fig. 1, from which it follows that both restored distributions practically 
coincide starting with a diameter of 10 um, gives an idea of the coordination of the data obtained in two ways. At the 
same time, one and a half dozen of the first points available for reading from a graph having a logarithmic size scale 
indicate the presence of a significant proportion of small, most dangerous for humans particles smaller than 10 microns. 
This circumstance is secondary since active air purification systems capture predominantly small particles. 

The fact, that at x => 10 um, the graphic distributions displayed by the Fritsch Analysette 22 Compact at 
different scales coincide, is confirmed by the following test. If we approximate both series of experimental data by the 
Rosin — Rammler curve and compare the obtained curves, we can quantitatively evaluate the effect of the digitization 


error on the final result. 
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Fig. 1. Digitization results of dispersion analysis data of metal-abrasive dust 


For approximation of the compared data presented in Fig. | by the model (2), the optimization problem related 
to mathematical programming was solved numerically [12]. The discrepancy between the actual and theoretical particle 
size distributions was minimized. If the discrepancy is determined by the Cartesian norm convenient for calculations, it 


is required to minimize the functional (®): 


(x, D,n) = Y[P(x,, D,n)-P.f > min, 
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where k is the number of the nodal (tabular) point; P, is the experimental probability density value; D and n are the 
described above parameters of the Rosin-Rammler model distribution. 

Solving this problem through Excel, we obtain the following result. For the data obtained from a graph linear 
in x, the values of D and n are equal to 48.37 um and 1.588 pm, respectively. Moreover, the curve approximating P(x, 
D, n) is characterized by an average residual with experimental points 0 00274 and a correlation coefficient of 0 894. 
Similarly, solving the problem for 35 right points reconstructed from a logarithmic graph, we obtain the values of D = 
48.61 um, 7 = 1.346, residual — 0.00187, correlation coefficient — 0.934. The proximity of the parameters D and n for 
both tables, together with a small residual and a high correlation, indicate the approximate equivalence of the studied 
graphic images of the analyzed result. The correlation of the functions P(x, D, n) with the calculated and above 
parameter values provides for the difference in the presentation of data in a formula form obtained from the compared 
graphic sources. Its value is 0.99. 

Since the dust dimensional characteristics averaged on the basis of a certain realistic distribution are important, 
we calculate the set of indicators Dij according to the Rosin-Rammler model initialized above for the values of the 
parameters D and n obtained in different ways. Then, we compare the results to each other and to the D, integral 
indicators recorded by the Fritsch Analysette 22 Compact. 

The purpose of comparing the analyzed graphic distributions to the integral indicators Dj recorded by the 
device is to find out exactly which distribution is shown on these graphs. The fact that on the screen of the Fritsch 
Analysette 22 Compact, the D4; value coincides with the arithmetic mean diameter in five decimal places, and the D3, 
value coincides with the geometric mean size just as accurately, suggests that not the parameter P(x,), but the 
differential x*-weighted distribution function is recorded: 

dF (x;) = P(xp) x0 dxp. (4) 
This is a key circumstance in interpreting the dimensional analysis results of dust particles based on histograms printed 
by the device. We verify the proposed hypothesis through calculating the moments Dj. The results of comparing the 
calculations when accepting or rejecting the hypothesis of the authors to the corresponding numerical data recorded by 


the Fritsch Analysette 22 Compact are given in Table 1. 






































Table 1 
Comparison of Sauter’s dust diameters based on approximations of digitized graphic data 
eecee Hypothesis (4) is not true Hypothesis (4) is true Instrument 

vit , D= 48.4 um, D= 48.6 um, D= 48.4 um, D= 48.6 um, numeric 

n= 1.588 n= 1.346 n= 1.588 n= 1.346 data, um 
Da 90.2 107.1 43.4 44.4 43.5 
Da 83.2 97.9 30.4 27.5 20.9 
Dy 75.2 87.1 15.6 12.5 8.18 
Dao 65.5 73.6 7.61 6.11 4.23 
Dsp 76.7 89.5 21.3 17.0 10.0 
D3 68.6 78.5 9.38 6.62 3.55 
Ds 58.9 64.9 4.26 3.15 1.95 
Dy 61.4 68.8 4.13 2.57 1.25 
Dy 51.6 55.2 1.90 1.35 0.66 
Dy 43.4 44.4 0.88 0.71 0.59 


























As can be seen from the table, the hypothesis on the essence of the data obtained using the Fritsch Analysette 
22 Compact service program is correct. This conclusion is crucial when interpreting the dust dimensional analysis 
results. 

Statistical analysis results of graphic data and their interpretation. To better describe the finely dispersed 
part of the dust, we use the graphical data (Fig. 1) in a logarithmic size scale and approximate them by a five-parameter 


distribution corresponding to the weighted sum of two lognormal distributions: 
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where D and D, are the peak position; 0), 62 are the peak width, o is the fraction of particles in the first mode. 

Solving the corresponding optimization problem using Excel tools yields the following result: D,; = 45.2 um, 
D, = 7.0 pum, 6; = 1.97 um, o2 = 2.1 um, a = 0.858. The average discrepancy between the approximating function and 
the initial data is 31%, and the correlation coefficient is 0.977. This is a satisfactory agreement given the high error of 
the experimental values (Fig. 1). Note that this result refers to the x°-weighted true size distribution function of dust. It 
can be interpreted as follows: the bulk of the dust (about 85%) is made up of particles larger than 10 um; therefore, for 
practical air purification, this distribution can be replaced by a two-parameter lognormal (2) one with the following 
parameters: D = D, = 45.2 um and o, = 06; = 1.97 um. However, such a simplification will not allow a qualitative 
approximation of D, by j < 3, while the consideration of the fine fraction is significant for some applications, for 
example, for calculating all the used diameters D, and the moments of the distribution function P(x). Below are the 
results of comparing the double lognormal and lognormal approximations of P(x). Fig. 2 shows the experimental data 


approximation by a double lognormal distribution: 








Igx-1.65 ies088 


aL . as 
P(x) = 0.176-e a 0416 + 0.563-e Al es 


The mean-square ratio error is 0.31, and the correlation coefficient with experimental points is 0.977. 
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Fig. 2. Approximation of experimental data by double lognormal distribution 


Fig. 3 shows the experimental data approximation by the lognormal distribution: 
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The mean-square ratio error is 1.23, the correlation coefficient with experimental points is 0.970. 
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Fig. 3. Approximation of experimental data by lognormal distribution 
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The possibility of approximating experimental data by the Rosin-Rammler distributions was also studied: 
e by classic generalized three-parameter: 
P(x, D, n, m) = (x/D)" - &!" / | (</Dy" - dx); 
e by simplified exponential: 
P(x) =1/D-e, n=1, 
Fig. 4 presents the experimental data approximation by the Rosin - Rammler distribution: 
P(x) = 7.73-10? x0 exp[(x/48.89)' >]. 
The mean-square ratio error is 0.47, and the correlation coefficient with experimental points is 0.972. 
Fig. 5 shows the experimental data approximation by a three-parameter distribution of the Rosin-Rammler 
type: 
P(x) = 6.66-10°-x°*-exp[(x/66.18)!7°]. 


The mean-square ratio error is 0.51, the correlation coefficient is 0.972. 


Rosin — Rammler unnormalized 


















































140 
° 
420 
408 
on 
Se 88 ¢ Experiment 
oy 6 — Approximation 
aw 48 
20 
-1 0 1 2 3 
Lg X, um 


Fig. 4. Approximation of experimental data by the Rosin — Rammler distribution 
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Fig. 5. Approximation of experimental data by three-parameter distribution of the Rosin-Rammler type 


It follows from the analysis that for the problems other than cyclone calculations, it is convenient to 
approximate the experimental data of the authors by a one-parameter monotonically decreasing function: 
P(x) =1/D-e*”. Fig. 6 shows the approximation of experimental data by the exponential distribution 


P(x)=1/47.13-e*""?. The mean-square ratio error is 0.42, the correlation coefficient with experimental points is 0.948. 
q 
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Fig. 6. Approximation of experimental data by exponential distribution 


Table 2 shows the integrated indices of dust dispersion calculated on the basis of model distributions in 


comparison to the evaluated software of the Fritsch Analysette 22 Compact. Model distributions are initialized by 

































































graphic data. 
Table 2 
Integrated indices of dust dispersion 
Mathematical model used to evaluate 
Ss | Be | 8 
3 5 g 2 SI s 
Parameter : & 5 g 5 5 Actual value 
ob o 1 a 1 lo} 
ro) ra go g oe 
4 = ze 5 iS 
fo) Oo . jo} 
lal ae la 
Da3, wm 51.7 44.8 44.7 45.2 43.6 43.5 
Da, wm 41.2 23.0 27.3 28.1 21.4 20.9 
D4, wm 32.0 9.45 12.2 12.9 8.61 8.18 
Dgo, pM 25.0 4.87 5.99 6.32 4.40 4.24 
D3, wm 32.8 11.8 16.6 17.5 10.3 10.0 
D3), wm 25.2 4.34 6.4 6.92 3.79 3.55 
D309, wm 19.6 2.32 3.07 3.28 2.03 1.95 
Dy, wm 19.4 1.59 2.46 2.74 1.39 1.25 
Doo, mM 15.1 1.03 1.32 1.42 0.904 0.858 
Dio, wm 11.8 | 0.666 0.704 0.738 | 0.586 0.588 
Mode, pm | 52.1 45.2 48.9 48.4 43.9 43.4 
MSR o, pm | 35.0 37.2 33.6 33.8 39.9 35.7 
Asymmetry | 1.49 1.17 1.19 1.20 1.41 1.25 
Excess 2.28 1.49 1.46 1.50 1.66 1.66 














The data in Fig. 4-6 and Table 2 provide for the comparison of the corresponding results to each other and to 
the approximation data of the lognormal and double lognormal functions. These results indicate a close approximation 
quality achieved using the five-parameter double lognormal, three-parameter Rosin-Rammler type, the two-parameter 
classical Rosin-Rammler type, and the one-parameter exponential model distributions P(x). Meanwhile, the lognormal 
distribution does not correspond to the experimental data since it does not reflect the presence of a significant numerical 
fraction of very small particles in dust. At the same time, this distribution is best suited for the practical purpose of the 


study — optimization of the air purification system. 
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The approximation by one-parameter monotonically decreasing function P(x) = 1/ D - e ” demonstrates good 
agreement with the experiment for small-sized dust fractions, which, due to their representativeness, significantly affect 
the mean-square ratio approximation error. 

Conclusion. The experimental data analysis performed by the authors regarding the size distribution of 
particles formed under a rail cutting reduces to the following: 

1. The Fritsch Analysette 22 Compact is not optimal for studying the fractional composition of dust since it 
does not present measurement data in primary numerical form. The digitization of the graphs provided by this device is 
the major source of error in the interpretation of the relevant data. 

2.For a detailed analysis of the dust particle size distribution based on the available experimental data 
(available from the authors), the simplest exponential particle size distribution P(x) = 1/47.13- e *4”*is best suited. 
Based on this distribution, all integrated indices provided by the instrumental measuring complex can be correctly 
reproduced, along with the graphical data. 


3. The tasks of air purification require using a classic lognormal distribution: 





pai 
(x) = 0.563-e : 
whose parameters are calculated by the authors using mathematical programming methods. They made it possible to 
calculate the important for designing cyclones value D3, = 32.8 um, which is three times higher than the numerical 
value produced by the Fritsch Analysette 22 Compact. This means that the application of the numerical values provided 
by the named device in the cyclone calculation will cause the inability to remove almost all large particles from the 
working area. 

4. Due to the fact that the efficiency of the cyclone is affected not only by the value of D3, but also by the 
density of the dispersed material consisting of steel and abrasive, a precise analysis is required to ensure the efficiency 
of air cleaning. Its subject is a separate study on the dispersion of metal and abrasive dust resulting from the operation 


of a rail-cutting machine. 
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Numerical modeling and experimental estimates of structural member 


fatigue characteristics 
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Introduction. In the algorithm for predicting the resource of machine parts, models of external actions, fracture re- 
sistance, and temporal development of a particular type of damage to these units interact. The applied issues of the re- 
sistance of machine parts to fatigue failure are considered. The scientific research and regulatory materials are adapted 
to determine the characteristics of endurance to the specifics of structures and materials of the construction-and-road 
machinery and agricultural machinery. The work objective is to use the analysis of existing methods to develop recom- 
mendations for the calculated determination of the endurance parameters of structural members of the agricultural ma- 
chines. 

Materials and Methods. The initial data were scientific studies on the fracture mechanics of the engineering materials 
and structures, as well as standards for endurance characterization methods. The need to test methods for determining 
fatigue characteristics to use them in the road construction machinery and agricultural machinery projects follows from 
the specifics of their designs, operation conditions and industry-specific range of materials. Based on the analysis of 
existing methods, it is required to develop recommendations for the calculation of the endurance parameters of structur- 
al parts of the agricultural machines. For this, the components of the load-bearing systems of the staged design were 
presented in the form of a set of plates of the corresponding thickness; and the concept of the critical radius of the stress 
raiser at the welding sites was also used. Numerical methods using mathematical models were applied. The calculation 
results were verified through comparing them to experimentally determined fatigue characteristics of the members of a 
combine harvester on a test bench. 

Results. Radius critical values of the stress raisers for various types of welds are obtained. The tables of the calculated 
and experimental endurance limits of the combine bearing components are well correlated and can be used under de- 
signing. 

Discussion and Conclusions. The theoretical foundations laid down in the study open up great opportunities for applica- 
tions to the design of various machines. The considered fragment of adaptation of the theoretical approach to the agri- 
cultural engineering objects can be used in the design of load-bearing systems in the related fields of engineering. 


Keywords: endurance range, fatigue fracture, fatigue crack, weld joint, computational and experimental methods. 


For citation: Yu.P. Man’shin, E.Yu. Man’shina. Numerical modeling and experimental estimates of structural member 
fatigue characteristics. Vestnik of DSTU, 2020, vol. 20, no. 1, pp. 79-86. https://doi.org/10.23947/1992-5980-2020-20- 


1-79-86 


Introduction. In the computational and computational-experimental methods for predicting the durability of 
the load-bearing units of ground-based mobile mass-produced vehicles (automobiles, tractors, road construction, irriga- 
tion and drainage, agricultural, military equipment, etc.), it is recommended to take the early formation of a fatigue 
crack of 0.2 — 0.5mm length under conditions of high-cycle fracture as the key reliability criterion’. In this case, the 
initial information on the structural properties is the characteristic of the fatigue curve. Consider the possibility of de- 


veloping methods for determining the fatigue characteristics of machine parts’. 


' Kogaev VP, Makhutov NA, Gusenkov AP. Raschety detalei mashin i konstruktsii na prochnost' i dolgovechnost': Spravochnik [Strength and dura- 
bility design of machine parts and structures: Reference Guide]. Moscow: Mashinostroenie; 1985. 224 p. 

°RP 206-86. Raschety i ispytaniya na prochnost'. Metody opredeleniya kharakteristik soprotivleniya ustalosti detalei mashin s uchetom rasseyaniya: 
Metodicheskie rekomendatsii [Calculations and strength tests. Methods for determining of fatigue resistance characteristics of machine parts consider- 
ing scattering: Recommended practice] Moscow: VNIINmash; 1986. 50 p.. 
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Structural parts are represented by plates working in push-pull. The stress concentration is caused by a sharp 
change in the thickness of the plates with a certain radius of curvature p at the junctions, up to a sharp cut. The stress 
concentration associated with the formation of a weld on the plates under study is considered the subject of independent 
research and is taken into account as a result. 

The correct parameter selection to fatigue damage of machine parts and load-bearing systems is the basis for 
calculating resources and predicting reliability [1]. Consider the methodology and the calculation results of the charac- 


teristics of the fatigue resistance of the elements (zones) of the loadbearing structure based on the model shown in 


Fig. 1. Scheme for the calculated determination of fatigue characteristics 


The model is a fragment of a welded joint of thick-walled units. It is typical of many self-propelled and trailed 
vehicles of land transport, for example, a frame system of self-propelled combines. The initial data for the calculation 
are given in Table 1. 


















































Table 1 
Initial data for determining fatigue characteristics of joint 
assemblies shown in Fig. | 
Option 1 2 3 4 
Steel grade 09G2 09G2 20 20 
Working corrosion corrosion 
conditions 
o,, MPa 460 460 420 420 
o,, MPa 300 300 250 250 
o_,, MPa 232 232 213 213 
Vo 0.069 0.069 0.062 0.062 
Kreg 0.85 0.55 0.85 0.58 
K, 0.50 0.50 0.50 0.50 
K, 0.95 0.95 0.95 0.95 
t, =t,,mm 10 10 4 4 
p,mm 0.1-2 0.1-2 0.1-2 0.1-2 
p,mm 0.5 0.5 0.5 0.5 
V, 0.36 0.36 0.36 0.36 
/,mm 200 200 80 80 

















The median endurance limit is calculated by the expression’: 


' Kogaev VP, Makhutov NA, Gusenkov AP. Op. cit. P. 112-224. 
? GOST 25.504-82. Raschety i ispytaniya na prochnost'. Metody rascheta kharakteristik soprotivleniya ustalosti [GOST 25.504-82. Strength calcula- 
tion and testing. Methods of fatigue strength behaviour calculation]. Available at: http://docs.cntd.ru/document/1200012858 (accessed: 


11.12.2019). 
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where G_, is the median endurance limit determined on smooth laboratory samples of standard diameter d, = 7.5 mm; 


K isthe endurance limit reduction coefficient: 











Guia 


Gig 


where K, = is the effective stress concentration factor; K,, is the full-size influence coefficient; K,,, is the sur- 


face condition influence coefficient; K, is the influence coefficient of surface softening from metallurgical welding 
processes; K_, is the anisotropy factor of material properties. 

For the example considered, we will use the recommendations [2]. We accept K, =0.5, For normal operating 
conditions, we accept K,,, = 0.85. During operation or storage of the machine with the ingress of freshwater, moisture 
and the formation of corrosion, we accept for steel 09G2 on GOST 5521-93 K,, =0.55, for steel 20 on GOST 2591- 
2006 K,, =0.58. The initial data for calculating the endurance of the samples and their variation coefficients V,_, are 


given in Table 1. 
To determine the characteristic relations of parameters, the technique’ recommends the formula: 


K, a, 
K 








do 0.5¢ 
where a, is a theoretical coefficient that determines the degree of increase in stress in the zone of their concentration, it 
is calculated as per the local stress-strain state using finite element models; € is a distribution parameter. 


The similarity equation for fatigue failure is close in form to the Weibull distribution. For plate elements, it has 


the form’: 


Zz (B+1)(B+2) Ce 


where G; is the gradient of the first principal stress, it is determined by the formula: 


t, +t 
J =-2.31g(1- p) == ( 
G 


1 





, (1) 


= B 
Gi =23/p; (+) = 0.0152(B+1), 
Zz 


here, V, is sensitivity of metal to the stress factor: 
V,, =0.2—0.0001o, ; 
1 
=—— 0.64. 
P V. 


o 


We obtain the median value of the parameter C with the probability of destruction P=0.5 from the solution to the 


similarity equation (1). For the parts in question, the values of C are given in Table 2. 


"RP 206-86. Op. cit. 1986. 50 p. 
> RP 206-86. Op. cit. 1986. 50 p. 
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Table 2 
Design parameters of the structural fatigue resistance (Fig. 1) 
Option 1 2 3 4 
P,»> mm 0.5 0.5 0.5 0.5 
Oc 2.2 2.2 1.7 1.7 
Vg 0.19 0.19 0.22 0.22 
C 2.12 2.12 2.51 2.51 
V, 0.154 0.154 0.158 0.158 
B 5.85 5.58 5.69 5.69 
Gi: 1/mm 4.6 4.6 4.6 4.6 
K,/ Ki, 2.07 2.07 1.35 1.35 
K 4.74 6.09 3.22 4.57 
0 0.49 0.49 0.20 0.20 
Vs nan 0.05 0.05 0.06 0.06 
Vinge 0.21 0.21 0.24 0.24 
AK,, , MPa 9.53 9.53 9.53 9.53 

















The average stress gradient in the zone of radiusp , which is the place of concentration, is determined for the 
condition: 
p=(p,,) ¥ 0.5—0.0004(6, - 500). 
For the steels under consideration, the critical value of the radius of curvature in the connection of parts p,,= 0.5 mm 
was established. It corresponds to the formation of fatigue cracks not running further '. By this criterion, the value a, 


was selected. Fig. 2 shows the dependences of the theoretical stress concentration coefficient on the radius of curvature 
for various geometric design parameters. The calculated median endurance limits for the parts under consideration are 


presented in Table 3. 
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Fig. 2. Dependences of theoretical stress concentration coefficient on radius of curvature: | is t, =t,= 10 mm; 2 is t, =t,=4 mm 





; Kogaev VP, Makhutov NA, Gusenkov AP. Op. cit. P. 112-224. 
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Table 3 
Design parameters to define joint assembly characteristics based on 
the initial data of Table 1 shown in Fig. 1 


























Option 1 2 3 4 
6-lg 49 38 66 46 
(o—lg)P=0.05 , MPa 32 25 40 28 
(o—lg)P=0.95 , MPa 66 51 92 64 
m 4.54 3:52 6.36 4.48 
N,,10° cycles 2.9 2.6 2.9 2.7 
yo 0.11 0.11 0.10 0.10 
o=1,,.5. MPa 34 : - 7 























The variation coefficient of the fatigue limit is determined from the expression: 


V ne = We omax + ae + Voc ? (2) 


where V__,, 18 the variation coefficient of maximum stresses in the concentration zone; V__, is the variation coefficient 


omax 


of the fatigue limit of the sample; V_, is the variation coefficient of the theoretical coefficient in the concentration zone. 


The variation coefficient V,_, is set' and presented in Table | as source data. The variation coefficient V,,,,, is 


determined from the expression”: 
— = = 2 
1+0”° 
where @ is the relative fatigue failure similarity criterion: 
—— 
88.3G, 





The design values of the parameters are given in Table. 2. The coefficient of variation V,, is determined by the recom- 
mendations *: 


0a, 
op 


Vo= 


ao 











p 
te 
To find the derivative module, a linear approximation of the function a, = e(p/ F), is carried out, which is 
shown in Fig. 2 in the vicinity of the average value of the parameter p. The calculation results are given in Table. 2. It 
also contains the results of calculating the variation coefficient of the fatigue limit from the expression (2). 
We determine the values of the fatigue limit under the assumption of its normal distribution with a probability 
of destruction of 5% and 95% from the formula: 


(9_4,)p =9_,,0+ZV4,), 
where Z, is the fractile of the normal distribution. Table 3 shows the fatigue limit values with the indicated fracture 
probabilities. 


Next, we present the dependences for calculating the fatigue curve parameters: 
e indicator of the slope of the left branch of the fatigue curve’: 


m=[5+S2))K,; 
80 


e fatigue curve point abscissa: 


195 = (0.997-0.003m) , 
N, =10°6"2 


¢ coefficient of sensitivity to cycle asymmetry”: 
W, =0.02+2-10~%6,. 





' Kogaev VP, Makhutov NA, Gusenkov AP. Op. cit. P. 112-224. 
? RP 206-86. Op. cit. 1986. 50 p. 
> RP 206-86. Op. cit. 1986. 50 p. 
4 RP 206-86. Op. cit. 1986. 50 p. 
> Kogaev VP, Makhutov NA, Gusenkov AP. Op. cit. P. 112-224. 
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The calculation results of these characteristics for the components under consideration are given in Table 3. 
The studies have established’, that in the zones of stress concentration, there are critical radii of curvature in 


the range of values p,,= 0.1-0.6 mm [1]. Atp <p, , there are other patterns of similarity of fatigue failure compared to 
the area p > p,, . In this case, a crack-like ultimately sharp notch is considered, which is characteristic of the weld zones 


of the joinable units in various engineering industries [3, 4]. Two fatigue limits are recommended’ to determine: 


1. Fatigue limit o_,, atp >p,,,. It is determined by the criteria for the appearance of the first microscopic fa- 


tigue cracks and depends on the radius of curvature. 


2. Fatigue limit of final fracture o_,,,, atp<p,,. It does not depend on the radius of notch curvature and is 


determined from the formula: 
_ 1L,5AK,, 


O_h exp =~. 
V t+ Pp Kp 
where AK,, is the fatigue crack growth threshold; ¢ is the step depth in the zone of the voltage concentrator. 


The threshold for the development of fatigue cracks for mild steels, considering the asymmetry coefficient, can 


approximately be determined from the formula *: 
AK, ~ 6, 74V1-R . (3) 


The parameter AK, for the units under consideration is indicated in Table 2 at R=1. The values o calculated from 


—lgxp 
the expression (3) are given in Table 3. 

To assess the appropriateness of calculating the fatigue of structural elements presented in the form of plates 
with different transition radii in the connection zones up to an ultimately sharp notch, bench tests of field samples were 


carried out. The test results are presented in Table 4. 




















Table 4 
Experimental characteristics of fatigue resistance of harvesting machine metal structures 
te (c.,, ), MPa (o_,) (o_,) N,-10° 
Weldaone P =0.05 MPa P =0.95 MPa ie. expels 
; . by by by by 
min |mean| max} min} mean | max : 
mean min mean max 
flange and 
drive axle 3.2 | 3.7 | 43 | 31 37 43 4 4 70 81 0.54 3.7 
housing 
gearbox 
preckotand: | 94. |. 39: | 49 |.at| “aor | ae] 94 19 54 65 | 0.23 | 3.6 
drive axle 
housing 
mounting 
plateand |) 35) 41|4s]33} 38 | 43] 21 18 55 62 | 027 | 39 
drive axle 
housing 
frame racks 4.2 | 49 | 5.6 | 38 45 52 35 29 55 64 0.14 4.1 















































For bench testing, the highest accuracy group was selected with a relative error 0.1< ¢ < 0.2 , confidence prob- 
ability B = 0.8. The number of tested structures was in the range of 11-14 units. 





' Kogaev VP, Makhutov NA, Gusenkov AP. Op. cit. P. 112-224. 
> RP 206-86. Op. cit. 1986. 50 p. 
> RP 206-86. Op. cit. 1986. 50 p. 
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The data of Table 4 show large scatter of empirical estimates. Even if we exclude the weld zone of the flange 
and wheel axle beam from consideration, whose spread of values is caused by the instability of the process, in a given 


confidence interval, the value o_, 2 differs by 1.8—2.9 times. 


Conclusion. A comparison of the calculated (Table 3) and experimental (Table 4) characteristics of fatigue re- 
sistance in the zones of welded joints shows their good compossibility. The confidence intervals of the experimental 
values overlap similar intervals of the calculated values for the areas of the wheel drive axle. The average values for 
normal operating conditions differ by 1.3 times, and coincide under corrosion. For frame metal structures, the calculated 
values under normal conditions differ 1.5 times, and under corrosion conditions, they practically coincide. Good agree- 


ment with real values is shown by values o calculated by the threshold of development of fatigue cracks. 


—lgxp 

The materials presented provide for the conclusion that there is sufficient life reliability of the calculated char- 
acteristics of fatigue resistance based on the proposed method at the forecasting stage. The correct characteristics of 
fatigue resistance with a known or obtained from a field experiment distribution can serve as the foundation for calcu- 
lating the design life with the required probability of failure-free operation resulting from the technical specifications 
for the corresponding machine [5]. 

Thus, through analyzing the existing methods, recommendations have been developed on the calculated deter- 
mination of the endurance parameters of structural elements of agricultural machines. Under the engineering change of 
a project or a product prototype up to standard durability through structural and technological techniques, the character- 


istics of fatigue resistance are used in methods of approximate and functional resource estimates [6, 7]. 
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Adaptation of structures of steel laser hardening zones to friction conditions 
A. V. Brover 


Don State Technical University ( Rostov-on-Don, Russian Federation) 





Introduction. The structural state of the surface layers of engineering products manufactured through laser processing 
under various irradiation modes is considered. Structures with the highest possible stability with respect to external ac- 
tion under friction conditions, or with the possibility of optimal adjustment and additional hardening during operation 
by the formation of secondary structures, were implemented. Under the operating conditions, under the impact of me- 
chanical and thermal pulses, an express rearrangement of one structure to another occurs, which is stable at a higher 
level of load-speed conditions. Thus, the phenomenon of structural-energy adaptability is realized. The resulting adapt- 
able structures most efficiently dissipate the energy introduced into the tribosystem, and minimize the wear of friction 
pairs. 
Materials and Methods. We studied samples of P6MS steel using the following equipment: 

e Kvant-16 processing station with a power density of 100 MW/m’ for pulsed laser irradiation of samples; 

e Neophot-21 optical microscope for metal physical studies; 

e DRON-O0.5 diffractometer for phase composition identification. 
Results. It is established that the material of the samples after laser treatment is able to efficiently dissipate the energy 
supplied during friction through its transformations at various structural levels. As a result, it becomes possible to con- 
trol the surface strength and wear resistance of materials using the concept of structural adaptability of friction pairs, 
which extends their range of performance. Wear resistance of the irradiated steels is determined by both their initial 
hardness and the work-hardenability during friction. It is established that the thermal-strength loading of the steel sur- 
face irradiated layers during friction destabilizes austenite to the ya transformation, i.e., it contributes to its transfor- 
mation into deformation martensite. 
Discussion and Conclusions. In relation to the specific loading conditions, it is required to regulate the amount and de- 
gree of stability of the residual austenite in laser-hardened steels and alloys, which provides the necessary operational 


properties. 
Keywords: laser irradiation, machine building materials, wear resistance, structural adjustability, surface strength 
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Introduction. To increase the efficiency of hardening technologies including laser processing, it is necessary 
to determine the possibilities of targeted use of internal reserves of structural adaptability of products for various func- 
tional areas from steel and alloys under operating conditions [1—8]. It is required to analyze features of the structure- 
energy state of laser-irradiated materials in friction pairs and determine major directions for increasing the wear re- 
sistance of tribosystems. 

Under laser processing, several ways of transforming the structure were used to goal-oriented changes in the 


properties of materials [9-12]: 
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e increasing the structure dispersion under local plastic deformation as a result of dynamic polygonization in 
austenite, the formation of microdomains (fragments) of high density dislocations inherited under accelerated cooling, 
and as a result of phase hardening during polymorphic transformation; 

e formation of nanoscale extractions under the impact of plastic deformation under thermal exposure; 

e development of deformation in irradiated zones of martensite under external temperature-force loading. 

It should be taken into account that structural transformation, on the one hand, is a mechanism of strain hard- 
ening determined by an increase in the martensite volume, and, on the other hand, causes relaxation of microstresses 
and additional development of plastic deformation. These two factors act simultaneously and are competing [13-15]. In 
the case of the predominance of the first factor, high strength with satisfactory ductility is provided in steels. If the se- 
cond factor dominates, the steel ductility increases significantly while maintaining a sufficient level of tensile strength. 

With an optimal combination of factors, martensitic transformation under loading provides the best complex of 
mechanical properties of the material [16, 17]. 

Research Objective. Problems of designing adaptable alloy structures with specified operational characteris- 
tics using laser irradiation under optimal conditions were considered. The following ways were used to dissipate the 
energy pumped through laser processing and during the subsequent temperature-force tribological action: 

e transformation of part of the energy into heat; 

e dissipation due to the motion of defects in the crystal structure and plastic deformation; 

e dissipation under structural-phase transformations at different scale levels. 

The last structural factor is crucial for improving the operational properties of products if the structures formed 
under laser exposure are adaptable to operating conditions [18-19]. Adaptable martensitic-austenitic structures with a 
predetermined ratio of components that are purposefully formed under laser processing should efficiently dissipate en- 
ergy in the friction zones due to internal transformations, especially when peak loads are reached. This is possible under 
laser irradiation under optimal conditions and is a prerequisite for increasing the reliability and durability of products. 

The study objective was to obtain scientific knowledge on the possibility of implementation of structures in the 
laser processing zones of steels and alloys that adapt to temperature-force loading under friction and thereby increase 
the wear resistance of products. 

Materials and Methods. Metal physical studies of samples from P6MS instrument steel were carried out using 
a Neophot-21 optical microscope, a DRON-O.5 diffractometer, and a PMT-3 durometer. Laser irradiation was carried 
out on a Kvant-16 pulsed-action facility with a radiation power density of 80-150 MW/m°. Wear tests were carried out 
on the MI-1M installation according to the “disk — block” scheme at a load of 500 N with a linear sliding speed of 190 
m/min. Samples in the form of disks went through a full cycle of volumetric heat treatment — they were quenched and 
triple-tempered. Some disks were irradiated along a side surface 10 mm wide. Inserts made of steel ShKh15 having 
hardness HB 130-180 served as counterfaces. Disk wear before and after laser irradiation was determined by weight 
loss through periodical weighing on an analytical balance with an accuracy of 10* g. 

Research Results. Metallographic and durometric analyses of hardened samples have shown that after laser 
treatment, a hardened layer with a depth of 80-120 um is formed on the surface. Fig. 1 shows the wear test results that 
demonstrate a distinct advantage of steels after surface laser treatment. Moreover, due to the running-in of friction pairs 


and structural changes on contacting surfaces with an increase in test time, the advantage increases by 13 times. 
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Fig. 1. Wear curves of samples from R6M5 steel after quenching at 1280°C 
and triple tempering (1), laser quenching (2) 
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Consider the features of structural and phase transformations realized in the surface layers of laser-hardened 
samples under friction and eventually causing an increase in the wear resistance of the material. X-ray diffraction stud- 
ies of laser-irradiated samples of P6MS steel were carried out before and after wear tests, for 15 minutes and 150 
minutes. The results are presented in Fig. 2 in the form of fragments of X-ray patterns taken from the surfaces before 
and after friction. 
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Fig. 2. Fragments of X-ray diagrams of R6MS steel: | is after laser quenching; 2 is after laser quenching and 
wear tests for 15 min; 3 is after laser quenching and wear tests for 150 min. 

Prior to wear tests, the surface layers of steel after laser hardening have a two-phase austenitic-martensitic 
structure (Fig. 2, curve 1). Moreover, the intensity of austenite reflections, and, consequently, its amount in the struc- 
ture, decreases significantly with increasing the test time (Fig. 2, curves 2, 3). This is explained by the fact that the con- 
sidered conditions of temperature-force loading destabilize austenite to martensitic transformation due to the precipita- 
tion of doped carbides from it. As a result of the y-phase depletion by carbon, it acquires the ability to transform into 
martensite deformations with the properties similar to cooling martensite. This provides for technological and structural 
plasticity of the steel surface layers (micro TRIP effect). 

It should be noted that martensite obtained through laser hardening, under the effect of temperatures and forces 
in the friction zones, also transforms due to strain hardening during wear (Fig. 2, curves 2, 3). Confirmation is the in- 
crease in the width of the martensite reflections in X-ray diffraction patterns caused by an increase in the density of dis- 
locations and the refinement of the blocks of its fine structure. As a result, the hardness and performance of laser- 
irradiated steel is further enhanced. Thus, under friction, the structure of the irradiated surface layers adapts to the load- 
ing conditions in the friction pair, which leads to an increase in the hardness and wear resistance of steels after laser 
treatment. 

Through changing the irradiation regime, a different amount of martensite and austenite was achieved in laser- 
hardened zones, and a connection between the development of phase transformations in the process of self-organization 
of structures and the properties of steels was established. The results of metal physical studies have shown that laser 
processing of P6M5 steel should be carried out with a power density in the range of 80-120 MW/m”. In this case, the 
maximum possible hardening of steel is achieved at the level of 10—11.5 GPa, which causes an increase in the wear re- 
sistance of irradiated surfaces. An additional factor in increasing the wear resistance of steel under friction without lu- 
brication is the preservation of insoluble carbides of tungsten, vanadium, chromium up to 30% in the laser processing 
zones. 

It is established that under the action of high temperatures and pressures in the contact zones of the laser- 
irradiated metal, additional precipitation of carbides of the strengthening action occurs [20, 21]. This is facilitated by the 
high density of dislocations in steels both immediately after laser processing and additionally arising during friction. 

Fig. 3 a shows the microstructure of irradiated P6MS steel after the wear tests. A reference line is plotted; 
along it, a histogram of the surface profile height distribution, shown in Fig. 3 b, was determined. The line crosses the 


dispersed precipitated carbides. From the histogram, it follows that carbides have sizes in the range of 2-10 nm. 
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a) b) 
Fig. 3. Structure of laser-quenched R6MS steel: a) optical microscopy results; 
b) histogram of surface profile height distribution 


The precipitation of the dispersed carbide phase oriented along the sliding lines as a result of friction tests is 
noteworthy. The texture of dispersed precipitates thus formed causes a decrease in friction losses under the maintenance 
of products for various functional areas after laser treatment. 

Discussion and Conclusions. We note the features of two-phase austenitic-martensitic zones of irradiated 
steels for different types of wear. The results considered in this paper were obtained under friction without lubrication. 
Under the abrasive wear, non-uniform structures of laser hardening have lower wear resistance compared to fully mar- 
tensitic structures, which is explained by intensive wear of the austenitic structural component. Under the conditions of 
adhesive wear, as well as during the lubricated friction, two-phase structures of laser hardening zones are characterized 
by high wear resistance. In this case, the martensite sections play the role of a strong frame, and the wear-out austenite 
sections form microcavities on the working surface that hold the lubricant and improve the friction conditions. 

Thus, to increase the wear resistance, it is required to design a structure adaptable to operating conditions with 
a certain ratio of the volumes of martensite, austenite and carbides by selecting the optimal laser processing mode. 
There are promising options for laser modification of the surface layers of steels with carbides or nitrides of doping el- 


ements, or laser quenching and subsequent tempering, the temperature of which depends on the steel grade. 


Conclusions: 
1. The possibilities of the purposeful use of internal reserves of structural adaptability of steels and alloys after 


laser surface treatment to improve their wear resistance under friction are determined. 

2. The analysis of structural-phase transformations in the laser-irradiated zones of steels under temperature- 
force loading of friction pairs has been performed. It is shown that the transformation of residual laser hardened austen- 
ite into strain martensite, strain hardening of martensite obtained by laser heat treatment, and the multiple precipitations 
of dispersed reinforcing carbides in the irradiated zones of steels, contribute to increasing the wear resistance of materi- 
als. 

3. Through changing the modes of laser treatment of steels, it is possible to form adaptable martensitic- 
austenitic structures in the irradiated zones with a given ratio of components that efficiently dissipate energy in the fric- 


tion zones due to structural transformations at different levels. 
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Introduction. The problem of automatic license plate recognition is considered. Its solution has many potential applications 
from safety to traffic control. The work objective was to develop an intelligent recognition system based on the application 
of deep learning algorithms, such as convolution neural networks that consider automotive standards for license plates in 
various countries and continents, and are tolerant to camera locations and quality of input images, as well as to changing 
lighting, weather conditions, and license plate deformations. 

Materials and Methods. An integrated approach for the problem solution based on the application of convolution neural 
network composition is proposed. An experimental analysis of neural network models trained to meet the requirements of 
the universal license plate recognition task was conducted. Based on it, models that showed the best ratio of quality and 
speed were selected. Quality of the system is provided through the optimization of various models with different 
modifications. In particular, convolution neural networks were trained using images from several datasets. In addition, to 
obtain the best results, the models used were pre-trained on a specially generated synthetic dataset. 

Results. The paper presents numerical experiments, the results of which imply the superiority of the developed algorithm 
over the commercial OpenALPR package on public datasets. In particular, on the 2017-IWT4S-HDR_LP-dataset, license 
plate recognition accuracy was 94 percent, and on the Application-Oriented License Plate dataset, 86 percent. 

Discussion and Conclusions. The resulting algorithm can be used to automatically detect and recognize license plates. The 
experiments show that the algorithm quality meets or exceeds quality of the commercial OpenALPR package. The 
developed algorithm quality can be improved through increasing the training dataset, which does not require the 
participation of the developer. 


Keywords: object detection and recognition, convolution neural networks, data generation and augmentation, license plate 
recognition 
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Introduction. Automatic License Plate Recognition (ALPR) systems are used for automatic speed control, 
identification of stolen vehicles, access control of vehicles in private premises and toll collection’ [1,2,3]. However, most 
existing algorithms” * work only for a specific license plate template or with complex image capture systems and are 
demanding on lighting conditions and types of vehicles [4,5,6]. 

Due to the rapid development of deep learning and its applications in the field of computer vision [7], it became 


possible to create an ALPR system capable of recognizing numerous license plate patterns in an arbitrary environment**' 





' Lotufo RA, Morgan AD, Johnson AS. Automatic number-plate recognition. IEEE Colloquium on Image Analysis for Transport Applications. Feb, 1990. 
P. 1-6. 

> Du S, et al. Automatic license plate recognition (ALPR): A state-of-the-art review. IEEE Transactions on Circuits and Systems for Video Technology. 
2013;23(2):311-325. 

> Gou C., et al. Vehicle license plate recognition based on extremal regions and restricted Boltzmann machines. IEEE Transactions on Intelligent 
Transportation Systems. 2016;17(4):1096—1107. 

4 Hsu G., Chen J., Chung Y. Application-oriented license plate recognition. IEEE Trans. Veh. Technol. 2013;62(2):552-561. 

> Li H., Shen C. Reading car license plates using deep convolutional neural networks and LSTMs. 2016. arXiv preprint arXiv:1601.05610. 
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[8]. The fundamental obstacles in recognizing license plates are specific shooting conditions — the presence of rain, snow 
or poor lighting. The recognition task becomes more difficult if the license plate has an uncommon area and aspect ratio, 
background color, shape, number of lines, font size, distance between characters, etc. 

The study objective was to develop a license plate recognition system that supports various regional license plate 
standards, and does not depend on the conditions of video recording of cars, such as dirt on the numbers, weather 
conditions, etc. 

Materials and Methods 

Automatic system of detection and recognition of license plates. In the framework of this section, we give a 
description of the software implementation of the system of detecting a vehicle in the photograph, identifying the license 
plate in it, and its recognition. 

The proposed system consists of a composition of several neural network models. The idea of the composition is 
that the result of the previous model is fed to the input of the next one, and the whole analysis process is divided into stages. 


Schematically, the principle of the system is shown in Fig. 1: 


Vehicle detection Licensee plane detection License plate recognition 





Fig. 1. Three-step approach for license plate detection and recognition 


The first stage of the algorithm is to detect a vehicle. To solve this problem, we used the SSD Resnet” model, 
trained on a COCO (Common Objects in Context) * dataset. The result of this model is fragments of the input image which 
include vehicles. At the second stage, the obtained fragments are fed to the input of the developed model for detecting car 
numbers, which returns image fragments that include only car numbers. At the final stage, the characters that make up the 
car number are recognized and glued together. A detailed description of the developed models used in the second and third 
stages of the algorithm is given in the respective sections of the paper. 

The selected approach has the following advantages: 

e objects that can be perceived as a license plate are cut off: signs in showcases, windows or fences; 

e it becomes possible to establish a relationship between the number and the corresponding vehicle providing for 
the system expansion in the direction of recognition of other characteristics, such as the type, make and model of the 
vehicle, color and direction of movement, tracking in the video stream; 

e system modularity consisting of independent models specializing in solving a specific problem at a high level is 


maintained. 


License plate detection model. The initial data set for training the model was 1700 images of Indian cars taken at 


a certain angle. The model trained on this set was not variable for different conditions; therefore, the pseudo-labeling* 





' Yuan Y., et al. A robust and efficient approach to license plate detection. IEEE Transactions on Image Processing. 2017;26(3):1102—1114. 

°Tensorflow detection model ZOO. Available at: http://download.tensorflow.org/models/object_ 
detection/ssd_resnet50_vl_fpn_shared_box_predictor_640x640_cocol4 sync_2018 07_03.tar.gz (accessed: 03.11.2019). 

3 COCO — Common Objects in Context. Available at: http://cocodataset.org/#home (accessed: 03.11.2019). 

“ Pseudo-Labeling and Confirmation Bias in Deep Semi-Supervised Learning / E. Arazo [et al.]. 2019. 8 Aug. // arXiv preprint arXiv:1908.02983. 
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method was applied for unlabeled data sets (images with vehicles). Thus, the initial data set for training was increased to 
400 thousand images. The size of the validation dataset is 5% of the training dataset. Images were randomly selected and 
checked for correct marking of license plate coordinates. 

To assess the quality of the model, IoU (Intersection over Untion), mAP (mean Average Precision) and AR 
(Average Recall) ' metrics were used: 

"‘DetectionBoxes_Precision/mAP': class accuracy averaged over IoU threshold values in the range from 0.5 to 0.95 
in 0.05 increments; 

'DetectionBoxes_Precision /mAP @.50IOU': average accuracy of classes by IoU value equal to 0.5; 

"‘DetectionBoxes_Precision/mAP (small)': average accuracy of classes for small objects (area <32’2 pixels); 

'DetectionBoxes_Precision/mAP (medium)': average accuracy of classes for objects (322 pixels < area <96%2 
pixels); 

"‘DetectionBoxes_Precision/mAP (large)': average class accuracy for large objects (962 pixels < area < 10000’2 
pixels). 

Table 1 


mAP result on validation data 


mAP 
| tak | Moget | mar Jrsrasoioy mAP@7510U fmAP (smath| (AP JenaP cng) 





Picense iat: || 8 MobileNen 69569 0.9843 0.9739 0.7199 0.8197 0.8544 
detection vl FPN 


The resulting model is resistant to the size of the input image (Table 1), detects license plates at an angle, as well as 
with an uncommon aspect ratio. 

Vehicle plate numbers recognition model. The problem of license plate recognition is reduced to the task of 
detecting objects of 36 classes: 26 letters of the uppercase Latin alphabet and 10 digits. 

The initial data set for training the model was 1700 images of Indian cars. Validation took place at 10% of the 
initial set consisting mainly of hard-to-read images marked up manually. 

For lack of a sufficient number of labeled data sets that meet the requirements of training the numbers recognition 
model, it was decided to generate a synthetic data set. Such transformations as rotation, blurring and darkening, as well as 


various sizes and colorings of the license plate, font styles, and the distance between characters were used as augmentations 
(Fig. 2). 


a 
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Fig. 2. Samples from generated datasets 





' mAP (mean Average Precision) for Object Detection. Available at: https://medium.com/@jonathan_hui/map-mean-average-precision-for-object- 


detection-45c121a31173 (accessed: 30.10.2019). 
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The model was trained on a synthetic dataset including 300 thousand images. To increase the accuracy and amount 
of real data for training the model, the pseudo-labeling technique was used, which increased the set of real data to 250 
thousand. Recognition accuracy on the validation dataset has increased by 15% reaching the mark of 63% (Fig. 4). This 


growth is due to the presence in real images of artifacts that are not available under generation. 


DetectionBoxes_Precision 


map 
tag: DetectionBoxes_Precision/mAP 


Fig. 3. mAP results of the model before and after adding real datasets 

















Table 2 
Comparison of mAP models on validation data 
Task Model mAP | mAP@,50IOU |mAP@,75IOU | mAP (small) mAP (large) 
as ode , ; sma fenedtin) arge 
Faster R-CNN | 9.4783 0.9372 0.3578 0.4689 0.5251 0.5543 
Resnet-101 
License 0.6332 0.9503 0.7228 0.6009 0.6397 0.7051 
plate 
va 0.5567 0.9534 0.6113 0.5569 0.6094 0.6493 
recognition | raster R-CNN 
Inception v2 
0.5784 0.9652 0.6184 0.5583 0.6157 0.6676 
































Research Results 
In the resulting combination of neural network models, each of the classifiers was independently trained using 


images from several datasets including license plates of countries from all continents and methods of increasing these sets to 
achieve stability under various conditions". 

IoU metrics were used to evaluate license plate detection results. For the task of recognizing license plate 
characters, the image in which all the characters in the license plate were correctly recognized was considered correctly 
marked out. An odd, missing or incorrectly recognized symbol was considered an error on the whole license plate. 

The results were compared to the latest cloud version of the commercial OpenALPR® package and the 2016 study* 
mentioned above. The testing was performed on the 2017-IWT4S-HDR_LP-dataset (Table 3) provided in the study” and the 


' Pseudo-Labeling and Confirmation Bias in Deep Semi-Supervised Learning. 

? Nowruzi FE, et al. How much real data do we actually need: Analyzing object detection performance using synthetic and real data. 2019. arXiv preprint 
arXiv:1907.07061. 

> OpenALPR Cloud API. Available at: https://api.openalpr.com/v2/ (accessed: 02.11.2019). 

‘Li H., Shen C. Op. cit. 

> Hsu G., Chen J., Chung Y. Op. cit. 
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Application-Oriented License Plate (AOLP) (Table 4) provided in the study', which includes 2049 images with Taiwanese 


license plates. They were divided into three subsets with different levels of complexity and shooting conditions: Access 


Control (AC), Law Enforcement (LE) and Road Patrol (RP). 







































































Table 3 
Numerical experiment on public 2017-IWT4S-HDR_LP-dataset 
Metrics Ours OpenALPR 
Number of correct license plate recognitions 619 / 653 377 / 653 
Percent of correct license plate recognitions 94.79 57.73 
Table 4 
Numerical experiment on public dataset 
Metrics Applied solutions AC Subset (%) LE Subset (%) RP Subset (%) 

ours 95.58 93.97 94.29 
OpenALPR 91.80 86.89 90.84 
[12] 1st approach(with CNN I) 93.53 89.83 86.58 
[12] Ist approach(with CNN I) 93.25 90.62 86.74 
IoU plate detection | [12] Ist approach(with CNN I 93.97 92.87 87.73 

& Ill) 
[12] 2nd approach(with global 90.50 91.15 83.98 

features only) 
[12] 2nd approach(with both 
local and global features) ee aes Beek 
ours 88.75 83.94 85.40 
Number of correctly 
recognized license 
plates: 

OpenALPR 86.04 77.98 85.71 
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Discussion and Conclusions. Models were built to detect and recognize vehicle license plates. According to the 
numerical experiments, the model for detecting numbers turned out to be more accurate than the solutions obtained in 2016 
and 2019" * on the dataset published in 2013°. The license plate recognition model turned out to be more accurate than the 
decision of 2019* on the dataset obtained in 2017 [8]. To increase the accuracy of the models, methods of generating 
synthetic data and pseudo-labeling were used. The resulting system has appeared resistant to the size of the input image and 


the environmental properties. 
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Introduction. A mathematical model of the industrial safety monitoring system in mechanical engineering is investigat- 
ed. The work objective was to create a mathematical model based on expert assessments of workplace safety parameters 
with a calculated and experimental justification of its applicability to the “STRAZH” expert security monitoring system. 
Materials and Methods. The classification of expert systems for engineering enterprises is proposed. The stages of cre- 
ating expert systems are considered. A methodology for assessing the consistency of experts as a basis for models of 
expert systems in the field of mechanical facilities safety is presented. 

Results. The basic safety parameters of the workplace are identified. A matrix of expert evaluation of parameters based 
on the opinion of leading experts in the field of engineering is created. The results of modeling the expert system 
“STRAZH” with the calculated and empirical support of the mathematical model validity are presented. The advantages 
of implementing expert systems to increase the level of personnel safety are proved. 

Discussion and Conclusions. The results obtained have a high degree of expert coordination and can be used in the de- 


velopment of expert safety monitoring systems for engineering enterprises. 


Keywords: modeling, expert systems, knowledge base, expert assessment, mechanical engineering, safety parameters, 


concordance. 
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Introduction. Difficult-to-formalize processes complicate full automation of engineering production. Smart 
manufacturing, smart enterprise and similar structures are implemented only on the basis of analysis information sys- 
tems (AIS) using the components of artificial intelligence. At present, a uniform, classical definition of smart manufac- 
turing is not presented in the technical and specialized literature. However, experts agree that “smart manufacturing”, 
“smart enterprise”, “smart factory” refers primarily to a widespread use of information technologies, computing devic- 
es, sensors and distributed networks to implement a highly efficient production process and provide its participants with 
maximum safety [1]. 

A modern approach to the development of intelligent AIS for smart manufacturing involves a widespread use 
of new methods for representing knowledge and programmed empirical algorithms for their processing [2]. 

First of all, among the AIS used in mechanical engineering, we single out two most promising classes. 

1. Management information systems (MIS) are designed for monitoring and management of difficult-to- 
formalize process facilities. Core components in the MIS structure include: 

— module for collection and processing databulk (big data) according to certain algorithms; 

— module for the expert assessment formation [3]. 

2. Expert systems (ES) are designed for the collection, processing and analysis of formalized experience of ex- 


perts in a specific field of engineering. Core components in the structure of ES include: 
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— module for the accumulation of expert knowledge in a specific field of engineering; 

— module for the formation of alternative control scenarios under specific conditions based on the empirical 
experience of experts [4]. 

AIS of both classes are complex software systems created to replicate empirical experience and algorithms de- 
veloped on its basis to increase the efficiency of engineering industries. 

The knowledge base is the central system component that is formed in the process of ES modeling, designing 
and operating. The main difference between ES and other information systems involves the solution to a clearly limited 
range of problems in a specific area [5]. Unlike traditional machine solutions, ES use not a procedural analysis, but the 
processing of deductive reasoning. Similar systems can find a solution to poorly defined and unstructured problems [6]. 

Materials and Methods 

ES in Mechanical Engineering. In the modern world, the accumulated, processed and analyzed knowledge is 
used for monitoring, preventing and forecasting emergencies that is the result of empirical studies of several generations 
of experts. In this regard, ES are essential in the modeling and prediction of dangerous events. 

MIS and ES are designed in two stages: 

— designing a module for the accumulation and structuring of knowledge in a specific field; 

— designing a module for developing recommendations and making a control decision based on specific facts 
and parameters for monitoring the state of an object. 

The use of ES in the field of labor protection at engineering enterprises is due to the need to reproduce the 
knowledge of experienced experts. This is one of the conceptual stages in the development of digital production. From 
the user point of view, ES are of interest at this time for a number of reasons: 

— they can solve various practical problems and in terms of results are not inferior to expert people; 

— they are focused on solving a wide range of tasks in nonformalized areas; 

— they do not require special programming skills, and working with them is available for a wide audience of 
qualified users [7]. 

In mechanical engineering, ES help to make decisions, manage facilities, identify emergencies and failures, 


and design production. Fig. 1 shows the basic classes of problems solved by ES in mechanical engineering [8]. 
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Fig. 1. Basic classes of problems solved by expert systems in engineering 


In the practice of machine-building industries for TV-7 machines equipped with a function for controlling the 
accuracy of product processing, ES of the “Archimedes 2008” type are used. Under processing, the base circles in the 
cross and longitudinal sections and geometric parameters are calculated using the “Archimedes 2008” system to identify 
possible deviations. At the same time, problems with deviation of the longitudinal section profile, deviation from 
roundness, ovality, are identified, errors of sizing, waviness, etc., are determined [9]. 

The experience of using ES in mechanical engineering made it possible to identify the main advantages of their 
implementation: 


- an increase in the quality of decisions made, 

- improving the quality of manufactured products, 
- increase in productivity, 

- advanced training of employees. 
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It should be noted that it is advisable to use ES to solve complex problems in engineering production [10]. 

The key concept of labor protection in mechanical engineering is “workplace”. This is the place where an em- 
ployee should be located or where he needs to arrive in connection with his job. It is directly or indirectly controlled by 
the employer. Workplace safety is regulated by the Occupational Safety Standards System SSBT (GOST 12). It should 
be noted that ES do not provide for full control of safety at the workplace. 

To increase the reliability of control decisions, a generalized expert assessment should be introduced into the 
workplace safety model. The key point in conducting an expert assessment is the selection of competent specialists with 
experience in the claimed field and capable of an adequate assessment of the technological situation [11]. 

Research objective is to develop a mathematical model and conduct a calculation and experimental justifica- 
tion of its applicability for ES security monitoring “STRAZH” (“System for the exact calculation of vital activity algo- 
rithms), based on the analysis of the subject area and expert assessments. 

Initial data. Based on the analysis of literary sources and the opinion of practitioners, 11 basic parameters of 
workplace safety were identified. 

. Equipment (functional content). 

. Compliance of the equipment with the anthropometric characteristics of an employee. 

. The availability of personal and collective protective equipment, as well as fire extinguishing equipment. 
. Access to the workplace and the ability to quickly evacuate. 


— 


2 
3 
4 

5. Serviceability of production equipment. 

6. Performing production operations in accordance with the requirements of technological documentation. 

7. Monitoring of distributed hazardous and harmful factors. 

8. Keeping the established order and organization, high production, technological and labor discipline. 

9. Qualification of the employee. 

10. Timely training and retraining of the employee. 

11. Regular monitoring. 

The totality of data on the key parameters of the workplace safety provides you for such a characteristic of the 
work process as labor intensity. This integrated characteristic of the labor process shows the load on the nervous sys- 
tem, sensory organs, and considers the emotional component. Labor intensity is normalized by types of loads: intellec- 
tual, sensory, emotional, monotonous, and operational. 

Development of mathematical model of the ES “STRAZH”. When developing the mathematical model of 
the ES “STRAZH”, 21 experts evaluated the safety parameters of the workplace on a scale from | to 12 points. The 
survey was conducted using questionnaires. Based on its results, a consolidated matrix for assessing workplace safety 
parameters has been created (Fig. 2). 
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Fig. 2. Workplace safety rating matrix 
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A key outcome of the peer review methodology is Kendall’s concordance coefficient, which measures the con- 
sistency of the expert group: 
W = 12 = (1) 
m -(n —n) 
where W is the concordance coefficient, m is the number of experts, 7 is the number of parameters, S is the sum of 
squared deviations of the rank sums obtained by each parameter from the average rank sum of ranks. 


The sum of squared rank deviations S is calculated from the formula: 
S= 3D? =3(d,-a);, (2) 
i=l i=l 


where D; is the rank deviation, i is the serial number of the parameter, d; is the parameter rank, d is the arithmetic mean 
of the parameter rank. 

The concordance coefficient varies in the range from 0 to 1: 0 corresponds to the complete inconsistency of 
experts, 1 corresponds to complete coordination. If the concordance coefficient is equal to zero, it is necessary to check 
the initial data and (or) analyze the membership of experts in order to replace them (partly or completely). If the coeffi- 
cient value exceeds 0.4—0.5, the quality of the assessment is considered satisfactory, if it reaches 0.7—0.8 — high. 

Thus, when calculating the concordance coefficient according to the formulas (1) and (2), we obtain the fol- 


lowing parameter values: 
S= v(d,-dy = 5476 + 2916 + ... + 361 = 39142, 
i=l 


12-S 12 -39142 


= = = 0.806. 
m(n—n) 22 -(1P-1)) 





Using Pearson’s “chi-square” criterion [12], the null hypothesis io: W=0 (expert opinions do not agree with 


each other), at an alternative h;: W #0 (expert opinions are consistent with each other) is tested. 


We introduce expert estimates, rank sums d;, rank sum deviations D; from the average d and D. in the design 















































Table 1. 
Table | 
The concordance coefficient calculation 
Experts 
Parameters} 1 | 2|}3]4/5]6]7]|8 | 9 |10/11)12]13)14]15]16)17|18/19/20/21 es ot Di 
1 11)10}11}11)/12}11)12)10)11) 8 }11}10/10}11)12)11)10)11)10)11) 7 221 74 5476 
2 10}11)12}10} 9 |10} 7) 9 | 7) 5 |10/11)12}10)10} 7 }12)10)11] 8 }10 201 54 2916 
3 9/9|]9}8}10)9}9 |11}9)9 |12)9]8 9|10}9}9/9]9/11 196 49 2401 
4 S/6/5)/7}/5/5/6/5)/5)/2]/5)5])6]}4/6)/5)/6)5)7]5 43 108 39 1521 
5 12)12}10/12}11)12)11)12)12)12) 9 }12}11)/12}11)12)11)12)12}12] 9 239 92 8464 
6 2/2)/2)2)4)2})2)3)2/]10)/2/2)2)2)3)2)2|2/2/2/6 58 —89 7921 
7 6/5/6}6}6/6/51)6 6/6/6}/5}6/5/6/3/6]/6)/6]5 118 —29 841 
8 414/4)4])/2)3|}4]4/4/4/4/4])4]5/4]4]/5/4/4/4/8 87 —60 3600 
9 7/7/8)5)7|8}10} 7 |10)7)7)7)7)7)7/9/8)7)5)7)4 151 4 16 
10 3/3/3}3)3/}4/3/2)3)3)3)3}3)3)2)/3)/4)3)3]3 412 72 —-75 5625 
11 8/8}7/9/817)8)8]8 11) 8} 8)9)}8)8) 8) 7] 8] 8 {10) 2 166 19 361 
1617 39142 
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S 4, 


The average rank sum of all parameters is ga2 se “ = 
n 





147. 


ee | 1 ; 
We use the expression d = a -‘m-(n+l)= - 21-(11+1) =147. as a control of calculations. 


To test the null hypothesis using Pearson’s “chi-square” criterion, we calculate the empirical value 
y° =m-(n—1)-W =21-10-0.806 =169.4, which is compared to the critical values of “chi-square” for the number of 


degrees of freedom n—1=10. 
The empirical value x” =169.4 falls into the critical region y7 > x5o,(n —1) (169.4 > 23.2), which allows us 


to reject the null hypothesis. The concordance coefficient differs significantly from zero; therefore, there is a fairly close 
consistency of expert opinions regarding the estimated parameters. 

Research Results. An ES is developed in three stages: modeling, design, construction [13]. At the modeling 
stage, an analysis of the subject area to identify the most significant links and relationships between objects is carried 
out; the totality of input and output parameters, the degree of their input on the processes under study are determined. 
To build a mathematical model of the ES “STRAZH”, the safety parameters of workplaces of the machine-building 
industries were identified. When assessing safety parameters, it became necessary to select empiric experts who were 
the most competent in the organization of the mechanical engineering processes, since there are no methods to guaran- 
tee single-value safety assessments. The experts selected were occupational safety engineers from leading enterprises of 
mechanical engineering in the Rostov Region, as well as leading lecturers from the Engineering Technology Depart- 
ment, Don State Technical University. 

Discussion and Conclusions. According to the study, the concordance coefficient reached 0.806. This indi- 
cates a high consistency of expert opinions. It is verified by Pearson’s criterion and is a prerequisite for the development 
of a high-precision ES model. 

In modern science, a significant place is occupied by the problem of decision support using ES. The introduc- 
tion of such systems in mechanical engineering will enable: 

— to reduce the time on solving complex security issues; 

— to reduce the likelihood of producing spurious solution; 

— to raise the level of labor safety. 

The study of this issue in the context of modern innovative production is of current interest. 
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Introduction. Fundamentally important problems of the structure mapping and transformation, functions and parameters 
of various properties in the systems of information knowledge objects organization, as well as functionally different pa- 
rameters in the management processes of transformations of material objects, are considered. For this purpose, structur- 
al-functional-parametric models and parametric models are respectively used. The distinction between these concepts is 
relevant and practically significant. The scientific novelty of the presented work involves studying methods and infor- 
mation tools used to determine functionally different technological schemes for the interaction of objects at the stages of 
design and pre-processing engineering. 

Materials and Methods. The concept of a “structural-functional-parametric model” is associated with the definition of 
the structure of the basic knowledge objects of the subject area. In this case, the “linking base” is the methods and cor- 
responding means of system engineering in the infological modeling technology that are used to solve practical prob- 
lems. The concept of a “parametric model” is associated with the solution to practical problems of the process control. 
Nature of these tasks is functionally different (technical, physical, chemical, biological). It should also be clarified that, 
in this case, we are talking about converting the parameters of various properties of real objects by methods and means 
of system engineering (almost a mathematical apparatus). 

Results. A “structural-functional-parametric model” and a “parametric model” are general theoretical concepts that have 
invariant properties necessary for solving practical problems of the subject knowledge area. Considering the organiza- 
tion system and management processes in this way, note that it is required to maintain data and logical connections be- 
tween them under static and dynamic settings. 

Discussion and Conclusions. To solve practical problems in the subject knowledge area according to the technology of 
information logical modeling, certain methods, tools, algorithms, and operations are used. The most complete mapping 
and transformation of information objects is possible only in structural-functional-parametric models and databases of 
their solutions. The application of structural-functional-parametric models is the most important condition for a success- 
ful transition to a high-level deterministic automation of information technology for solving practical problems of the 
subject area. As an example of such a problem, we can cite the machining production design engineering. 


Keywords: production design engineering, cutting, system analysis, information technology, decision modeling, system 
technique. 
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Introduction. The development of systems engineering allows us to distinguish between the concepts of “par- 
ametric model” and “structural-functional-parametric model” taking into account the properties of the concepts of 
“mathematical modeling” and “logical modeling”, as well as the functions of their knowledge [1, 2]. Systems engineer- 
ing technology for determining the structural-functional-parametric model is based on the formalization of domain 
knowledge. Formalized concepts are first built into the structure of the basic knowledge objects, and then into the struc- 
ture of basic knowledge objects of a higher level of the subject area. At this, the semantic and syntactic properties of the 
concepts [3] should be as strict as possible. In this case, the author’s updated interpretation of existing and new concepts 
is fundamentally important. By actualization, we mean activation, initialization of meaning by transferring from a static 


(outdated) state to a dynamic (updated) state while preserving all the necessary links and relationships within the system 
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and with the external environment. This approach to the formalization of knowledge of the subject area provides for its 
transfer to a higher level of organization. 

We introduce the updated concept of “system principles”. These are information logical statements of the ap- 
proach to solving the problems of researching manufacturing systems and processes of functionally different purposes. 
These principles are verified by the multiple practice of their use in the setting environment (static — for systems and 
dynamic — for processes). 

Materials and Methods. We formulate system principles. 

1. Information structural-functional-parametric and parametric models are investigated. The distinction is 
based on two system principles: 

- unity and generality, 

- difference. 

Let us compare the structural-functional-parametric models of information objects and the parametric models 
of real objects. 

Their unity and generality are manifested in the functions of destination under the computer processing of in- 
formation: display, transformation, storage, transmission. 

The difference is manifested in the methods and means of solving practical problems of the subject domain de- 
fined in the considered models. 

Technologies for determining structural-functional-parametric models are implemented by methods and sys- 
tems engineering tools: systems theory, systems analysis of decision making, set theory, and graph theory. 

Set theory uses graphic means: formalized notation of concepts, logical operators for making statements by 


imposing links between formalized concepts. 


Graph theory specifies overlays, intersections, and associations of links between functionally unified elements 
of the object structure according to the semantic and syntactic properties of the formalized concepts used. To this end, 
structure graphs and Venn diagrams are used. 

To determine the algorithms in the technology of modeling systems for organizing information objects, con- 
ceptual notions are used: mapping, transformation, structure, set, etc. 

In the practice of solving problems, systems engineering technologies are widespread. They are used in the 
processing of information on computers and servers (that is, in computer systems), as well as in the management of 
functionally different organization systems. 

In [4], systems engineering technology is presented as system technique technology. This significantly narrows 
the scope of its application. The paper [5] can be recommended for information on the methodology of systems engi- 
neering in information technology for identifying large, sophisticated powerful systems in solving practical problems. 

Technologies of system techniques are implemented by various mathematical methods and means. With their 
help, parametric models, based on the mapping and transformation of the parameters of various properties of real ob- 
jects, are determined. In the parametric models, processes of controlling the transformation of real objects within the 
corresponding domain organization systems are mathematically modeled. 

The authors of [6] investigate the solution to problems of the practice of manufacturing preparation of machin- 
ing production. Here is an example of significant difficulties that arise under the implementation of this task by modern 
means of the mathematical apparatus on an insufficiently formalized reference basis of knowledge. 

In [7], the practice of successful implementation of the systems engineering technology in functionally differ- 
ent management is described. In [8—13], various approaches to systems engineering technology for solving problems of 
controlling the mathematical apparatus are considered. The papers [14—15] providing for determining the state of the 
manufacturing preparation knowledge, are directly related to solving the problems of its practice. 

The purpose-oriented properties of the “structural-functional-parametric model” concept and the 
knowledge functions following from its formulation. In the further presentation of the material, we use the updated 
concepts introduced in [1]: 


- structure, 
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- main structure elements (integration, disintegration), 

- design quality of the main elements (integration, disintegration) in the structure of objects {components} of 
the DPP stage and the structure of objects {initial blanks} of the PPE (TO) stage, 

- basic object of knowledge, 

- information “slave” (main) transformation object (TO). 

The concept of “structure” will be considered in various aspects on the example of manufacturing preparation 
of machining production. 

We introduce the formalized notation of concepts: 

- DPP is the stage of designing preparation of production, 

- PME (TO) is the stage of preproduction machining engineering, 

- {...} 1s multitude. 

1. Consider the well-known judgment: “Everyone sees the material, the content is found only by someone who 
has something in common with it, and the form remains a mystery for most ... The form needs to be digested as well as 
the material, but it is much more difficult to digest”’. 

Between the elements of knowledge in the triad of cognition, the sequence of their link (material > content > 
form), since the content can change depending on the perception of the material, and for the same content the form of its 
display can change. 

In most cases, knowledge is developed through transforming their display form, which determines the level of 
achievement of their depth and changes accordingly. A necessary condition for the transition from the knowledge con- 
tent to the form of their display is the achievement of a certain level of combining concepts into a system. 

2. A single environment of one level of the structural elements of the “slave” transformation object, in a gen- 
eral case, interacts (is interfaced) with a group of structural elements of the “leading” object. This interaction is provid- 
ed by the imposition of functionally different relationships along the “reference” interface points from the side of the 
“leading” object structure to the structure elements of the coordinate systems of the “slave” object. It is on this coordi- 
nation that structural-functional-parametric models are based. 

3. The connection of the “slave” and “leading” objects of interaction enables to introduce the following con- 
cepts to define the concept of “structure” of basic knowledge objects: 

- composition of the structural elements, 

- type of working (functional) connection, 

- communication overlay method, 

- functions of working relationships, 

- parameters of the properties of structural elements and the relationships between them. 

4. The connection of objects is based on the system principle formulated in [16]: “The boundaries of their con- 
nection carry the greatest information about sets”. The structural elements of each of the interaction objects, determined 
by the boundaries of their connection, are functionally new. In this case, the tasks on determining the position of the el- 
ements in the coordinate system of the “slave” object are solved. This connection of objects is used to determine the 
structure of the basic knowledge objects of any subject area. 

5. The existence of the structure of the base knowledge objects is provided only by those initial (previously de- 
fined) and obtained new concepts that can be defined as formalized (symbolic properties of concepts [3]), as well as 
unified analogues of structural elements and the relationships between them, to be embedded in the structure. 

6. The main first and second base knowledge objects of the first type of the PME (TO) stage are specified in 
the corresponding two process charts for the interaction of “slave” objects according to the working functions per- 
formed by them (basing, geometric shaping of the form element). 

7. The main base knowledge object of the fourth type of DPP stage is determined by the type of the main first 
base knowledge object of the first type of PME (TO) stage in the corresponding process chart of the interaction of ob- 
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jects. In the structure of objects {components, assembly units} of the DPP stage, it is possible to identify three more 
“slave” objects according to the work functions performed by them (basing, guides, torque transmission, division and 
fixing) [17]. 

8. The scientific novelty of the presented paper involves studying methods and information tools used to de- 
termine functionally different process charts for the interaction of objects of the PME (TO) and DPP stages. 

9. The information technology for the automated solution to the PME (TO) tasks assumes that the most im- 
portant condition is met — the definition of two “slave” transformation objects as components in the structure of the 
{initial blanks, blanks} of the PME (TO) stage and in the general case of four “slave” transformation objects as an inte- 
gral part in the structure of objects {components, assembly units} of the DPP stage. 

10. The model information of the main base knowledge objects of PME (TO) is distributed in two parts: invar- 
iant and typical object-oriented parametric ones. They are logically unified and cannot be considered separately. 

The invariant parts of the information model structure of the main base knowledge objects at the PME (TO) 
stage are designed to solve the problems of optimizing material and labor resources in the main and auxiliary technolog- 
ical operations. Typical object-oriented parametric parts of these main base knowledge objects serve as the basis for de- 
termining the base knowledge objects of a higher level of the subject area of the PME (TO). 

11. Each of the two parts in the structure of information models of the main base knowledge objects is deter- 
mined by the form of display (graphic in two-dimensional space, the corresponding analytical in three-dimensional 
space). This provides the conditions for the automated conversion of the form of their display from one view to the cor- 


responding another and vice versa. 


12. The information content and form of the base knowledge object can be determined due to its structural- 
functional-parametric information model for displaying and transforming the algorithm into technology for solving 
practical problems of the subject knowledge area from input to output under computer processing. Such information is 
perfect. 

13. The main base knowledge objects (of the fourth type of the DPP stage, of the first and second types of the 
of the PME (TO) stage) form a systemic fundamental for the hierarchical structure of the knowledge base of the PME 
(TO) stage according to the seven levels of classification of its base objects of various types. Base knowledge objects of 
all types are defined in the contours of the coupling circuit of their structure elements according to the superimposition 
functions (relations, links) of the structure elements in the totality of logically and information-related coordinate sys- 
tems. The fundamental of the knowledge base structure of the PME (TO) stage is made up of elementary and composite 
superimposed simple and complex functionally different{elements of geometric shape}. These {shape elements} are 
used primarily to determine the elements in the structure of the group of objects of interaction of functionally different 
flowcharts, and secondly — to determine elements in the structure of the main base knowledge objects of the first and 
second types of the PME (TO) stage, as well as of the fourth type of the DPP stage. The first first-type base knowledge 
object of the stage is {basing flowcharts} on their possible set for fulfilling the functions of basing production objects 
and cutting tools into machining attachments, as well as the devices — into the corresponding working bodies of the 
cutting machine. The second first-type base knowledge object is {basing flowcharts of the shape elements into working 
machines and geometric formation shape elements by cutting on working machines} on their possible set. The second- 
type base knowledge object is {working machines} as an information and logical connection of the first first-type base 
knowledge object in a general case with a limited necessary set of the second first-type base knowledge objects. The 
third-type base knowledge object is {working machine systems} as an information and logical connection of the fourth- 
type base knowledge object of the DPP stage in a general case with a limited necessary set of logically and information- 
related second-type base knowledge objects of the PME (TO) stage. The fourth-type base knowledge object of the DPP 
stage is {basing flowcharts} on their possible set for fulfilling the basing functions for parts and assemblies into design 
products. The fifth-type base knowledge object of the DPP stage is {design products}. 

Research Results. In the paper presented, structural-functional-parametric models of functionally different 


types are used to determine the hierarchical structure of the knowledge base of the PME (TO) stage at seven levels. The 
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structure of the hierarchical knowledge base of the PME (TO) is determined from the results of extensive practice of 
cooperation with metalworking enterprises. This is the most complete structure of the composition of the elements and 
relationships between them. It has invariant properties with respect to a possible set of subject knowledge areas for its 
distribution, for example, pressure, welding, computer processing of information, management. 

Discussion and Conclusions. The hierarchical structure of the knowledge base of the PME (TO) subject area 
is grounded on: 

- structural-functional-parametric models of the main base knowledge objects of various types; 

- databases formed on a possible set of model solutions. 

Base knowledge objects of all types are distributed among the levels of the knowledge base structure. These 
objects are connected through an organic unity of the design quality parameters of the main elements (integration, disin- 
tegration) of their structure [1], which makes their separate consideration impossible. 

The solution to practical problems in the invariant parts of the structure of information models for the main 
base knowledge objects provides for an adequate calculation of the material and human resources needed to optimize 
the design and organization of the CKD production. 

The solution to practical problems in typical object-oriented parametric parts of the structure for information 
models of base knowledge objects of all types provides for the design quality of the basic elements (integration, disinte- 
gration) of the structure of “slave” transformation objects. 
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